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20.  >^bstract  (Continued) 

These  areas  of  sparse  data  in  the  WMO  compilation  have  been  covered  in  this 
investigation.  The  tentative,  first  estimate  of  the  distribution  of  lightning  dis- 
charges over  a large  area  of  Earth  differed  significantly  from  the  distribution 
of  thunderstorm  days  for  April  and  July,  but  the  two  distributions  were  simi- 
lar for  January  and  October.  Centers  of  relatively  high  occurrence  of 
lightning  discharge  on  a yearly  basis  were  located  over  South  Africa,  the 
Mediterranean  Sea,  Arabia,  Southeast  China,  Southeast  Asia,  and  Australia. 
The  occurrence  of  lightning  discharge  was  shown  to  follow  Sun  northward 
from  January  to  the  maximum  poleward  thrust  in  July.  The  occurrence  of 
lightning  discharge  receded  equatorward  from  July  to  the  end  of  fall  in 
December.  The  most  reliable  planetary-scale  estimate, of  the  average  inci- 
dence of  lightning  discharge  was^.  2 X lO"^  km'2  sec"^  ior  the_^rthern 
Hemisphere  (0-179E)  in  Januar^972  and. 3. 0 X 10_'5,km"“_,3ec"^or  April 
1972.^There  was  evidence  tha^  in  general,  the  areal  concentration  of 
lightnmg  discharge  in  the  Nor^nern  Hemisphere  (0-179E)  decrease  with  geo- 
graphia  latitude. 
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The  Distribution  of  Thunderstorm  Days, 
Lightning  Discharges,  and  the 
Incidence  of  Lightning  Discharge 
Derived  From  VLF  Sferics  Data 


1.  iMBoi)M;no> 

l.l  Nci-d  for  I his  Intrslifculion 

•Scientist.s  have  .studied  electrical  activity  in  the  atmosphere  for  more  than 
200  years,  ^ yet  there  remain  gaps  in  our  knowledge  of  the  natx  re  of  these  complex 
phenomena.  The  lightning  discharge  is  perhaps  the  most  dramatic  electrical  event 
in  nature  and  one  of  the  most  misunderstood.  ^ There  is  a need  to  understand 
better  all  aspects  of  the  lightning  discharge,  including  its  climatologv. 

3 

Horner  stated  that  there  is  a requirement  in  radio  science  for  a better  under- 
standing of  the  temporal  and  spatial  distribution  of  the  areal  density  of  lightning 
discharge.  He  wrote: 

If  a more  complete  knowledge  existed  of  the  densities  of  lightning 
discharges  in  different  parts  of  the  world,  it  could  be  coupled  with  a 
knowledge  of  the  energy  radiated  by  an  average  discharge  to  give  the 
world  distribution  of  radiated  power.  Hitherto  this  distribution  has 


(Heceived  for  publication  17  May  1977) 

1.  Chi-C'hen,  I,.,  and  Chi-Zhang,  I..  ( IDliG)  AFCK I.,  Contract  AF19(r>28)-5073, 

Emm-fi7-16.')  Translations,  Translated  from  Acta  Meteorologica  Sinica 
3G;275-279, 

2.  Viemeister,  P.  E.  (19fil)The  I.ightning  Book,  Doubleday,  New  York,  31fi  pp. 

3.  Horner,  F.  (1994)  Advances  in  Radio  Research,  Vol.  2,  Academic  Press, 

New  York,  pp.  lil -204. 
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been  assumed  to  follow  the  distribution  of  thunderstorm  days,  for  want 
of  a better  index. 

4 

Martin  and  Hildebrand  have  expressed  the  need  for  a climatology  of  thunderstorm 
and  lightning  parameters,  as  stated; 

Many  scientific  and  eng'ntering  disciplines  have  an  interest  in  the 
geographic  distribution  of  thunderstorm  activity  ...  In  much  of  this 
work  there  is  need  for  information  regarding  the  expected  thunderstorm 
activity  at  various  geographic  locations  around  the  world  to  enable 
assessment  of  its  contribution  to  the  total  atmospheric  noise  level  at 
any  other  geographic  location.  On  a worldwide  basis,  knowledge  con- 
cerning the  number  and  location  of  lightning  discharges,  which  create 
large  amounts  of  radio  frequency  energy,  is  still  very  meager  and  based 
to  a large  extent  upon  studies  made  in  1925. 

5 

Ryers  explained  that  meteorologists  have  a special  interest  in  atmospheric,  elec- 
trical phenomena: 

The  main  feature  that  distinguishes  meteorology  from  the  other 
sciences  that  might  concern  themselves  with  the  atmosphere  is  its 
emphasis  on  circulations.  These  circulations,  ranging  from  small 
turbulent  eddies  to  the  air  flow  over  the  planet  as  a whole,  influence 
the  distribution  of  electrical  properties.  Thus,  a two-way  relation 
between  lightning,  as  well  as  all  other  electrical  phenomena,  and  mete- 
orological conditions  exists. 

A better  knowledge  of  the  distribution  of  lightning  discharges  would  be  useful  in  the 
planning  of  aircraft  routes,  especially  over  remote  areas. 

Extensive  data  are  available  on  some  aspects  of  the  spatial  and  temporal  occur- 
rence of  thunderstorms.  There  is  a relative  paucity  of  information  on  the  occur- 
rence of  lightning. 

1.2  l.ili'rnluri'  Krvirw:  .Slatislic.s  of  rhundiT'.tomi  and  l,i;;h(ning  I’aranirlcrs 

A thunderstorm  day  is  defined  as  a local  day  on  which  thunder  is  heard  at  a 
weather  station.  The  observation  of  lightning  without  the  sonic  noise  that  it  pro- 
duces (thunder)  is  not  a sufficient  criterion  for  a thunderstorm  day  to  be  recorded.^ 
The  thunderstorm  day  is  the  only  planetary-scale  estimate  of  thunderstorm  occur- 
rence that  is  available.  Although  a useful  statistic,  it  has  obvious  deficiencies. 

4.  Martin,  J.N.,  and  Hildebrand,  V.  W.  ( 19(>5)  NOl.C  Kept.  628  Project  No. 

WK5-0062,  Naval  Ordnance  Laboratory,  Corona,  Calif.,  95  pp. 

5.  Hyers,  H.  R.  (1965)  Problems  of  Atmospheric  and  .Space  Electricity,  Elsevier 

Pub.  Co.,  Amsterdam,  pp.  491-496. 

6.  World  Meteorological  Organization  (1953)  Part  1,  Publ.  No.  21,  TP  6. 


Thunderstorm-day  data  provide  no  estimate  of  diurnal  variation,  duration,  or 

electrical  intensity  of  thunderstorms.  A thunderstorm  day  is  recorded  at  a station 

whether  there  occurred  a single,  isolated  thunderstorm  or  several  thunderstorms 

during  the  day.  Moreover,  the  efficacy  of  thunderstorm-day  statistics  is  based  on 

the  obviously  false  premise  that  thunder  always  is  heard  no  matter  where  it  occurs. 

Personnel  of  the  World  Meteorological  Organization  (WMO)  compiled  the 

thunderstorm-day  statistics  from  raw  material  submitted  by  the  meteorological 

services  of  the  member  nations.  The  following  quotation  is  taken  from  Part  2 of 
7 

the  data  presentation. 

It  should  be  made  clear  that  this  project  is  not  a detailed  climato- 
logical study  of  thunderstorm  activity  over  the  world  . . . The  maps  can- 
not be  considered  to  be  in  any  way  final;  they  are  subject  to  revision  in 
light  of  new  data. 

Brooks,  in  a classic  paper,  pointed  out  that  thunderstorm-day  records  may 
be  inaccurate.  He  said: 

Thunderstorms  which  pass  directly  over  the  station  may  be  noted, 
but  those  which  occur  at  a distance  of  several  miles  are  often  ignored; 
this  is  especially  the  case  in  tropical  stations  where  thunderstorms  are 
severe  but  extremely  local  — at  certain  times  of  the  day  in  the  rainy 
season  distant  thunder  is  so  common  that  it  simply  does  not  occur  to 
the  observer  to  enter  it  in  the  register  — in  fact,  he  may  not  be  con- 
sciously aware  of  its  occurrence  . . . 

In  the  same  paper.  Brooks^*  estimated  from  a review  of  station  reports  that 
there  arc  approximately  1800  thunderstorms  in  existence  on  Karih  at  any  one  time 
and  that  roughly  100  discharges  occur  per  second.  Heydt  and  Volland^  used  a 
heterodyne  type  of  receiver  to  study  the  amplitude  spectra  of  received  atmospherics. 
They  estimated,  based  on  counts  of  signals  at  frequencies  of  5 kHz,  10  kHz,  and 
40  kHz  that  approximately  120  discharges  occur  per  second.  This  result  is  in 
close  agreement  with  that  of  Brooks,  especially  when  one  considers  the  disparate 
methods  of  analysis. 

Aiya**^  used  lightning  flash  counters  to  estimate  that  a local  thunderstorm  in 
India  lasts  3 hours  on  the  average  and  produces  about  one  discharge  per  km^.  This 
corresponds  to  an  incidence  of  discharge  of  approximately  9 X 10  km  sec  . 


7.  World  Meteorological  Organization  (195fi)  Bart  2,  I^ubl.  No.  21,  Tf^  li. 

8.  Brooks,  C.  K,  B.  (1925)  Met.  Office  Geophys.  Mem.  and  Brofl.  Notes,  No.  24, 

Isjndon,  147-154.  ' ~ 

9.  Heydt,  G. , and  Volland,  11.  (1954)  J.  Atmos.  Terr.  Bhys.  25:85-104. 

10.  Alya,  S.V.C.  (1958)  Electro-Technology,  J.  .Soc.  Electronic  Engineers 
(Bangalore)  12:(No.  1):1-12. 
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Hornpr^^  reviewed  data  based  on  various  sources  (for  example,  lightning  flash 

counters,  strikes  to  power  lines,  and  radio  noise)  and  estimated  that  the  incidence 

- 5 - 2 - 1 

for  the  "main  thunderstorm  areas"  of  the  world  is  10  ‘ km  sec  . He  also 
estimated  that  the  incidence  for  the  world  as  a whole  is  of  the  order  10  ^ km  ^ 
sec"^.  Neither  Aiya  nor  Horner  indicated  over  what  period  the  respective  esti- 
mates of  incidence  were  averaged.  By  assuming  that  on  the  average  one  thunder- 

storm  occurs  per  thunderstorm  day,  one  may  infer  from  Brook's  work  that  the 

-4  -2  -1 

incidence  of  discharge  is  of  the  order  of  10  km  sec  for  the  world  as  a whole. 
1.  2.  1 THE  l.IGHTNING  DISCHARGE 

I.ightning,  in  general,  is  the  visible  electrical  discharge  produced  by  thunder- 
storms. The  lightning  discharge  is  a series  of  electrical  processes  by  which 
charge  is  transferred  between  centers  of  opposite  polarity.  I.ightning  discharges 
almost  always  are  accompanied  by  thunder.  In  April  1855,  five  "brilliant  flashes" 

of  lightning  were  seen  to  strike  the  Washington  Monument  without  any  observation 

2 

of  thunder.  However,  observations  of  this  nature  are  extremely  rare.  The  types 
of  discharge  which  are  within  the  scope  of  this  study  include  cloud-to-ground  and 

cloud  discharges.  Other  types  of  discharge  (for  example,  cloud-to-air,  and 

13  2 

cloud-to-cloud)  occur  very  rarely.  ' ’ 

The  cloud -to-ground  discharge  is  a composite  and  complicated  event.  It  may 

be  studied  conveniently  in  three  stages:  (1)  The  initial  stage  of  the  step  leaders 

forms  the  conductive  path  to  ground  and  ends  at  the  first  return  stroke;  (2)  the 

intermediate  stage  includes  all  return  strokes;  and  (3)  the  final  stage  comprises 

14 

the  residual  variations  in  electric  field  after  the  last  return  stroke. 

The  leaders  follow  a stepped,  tortuous  trajectory  from  the  cloud  to  the  ground 

and  ionize  a conductive  path  which  typically  connects  the  negative  center  of  charge 

in  the  cloud  to  an  induced,  positive  center  on  the  ground.  The  return  stroke  gives 

off  light  as  it  apparently  moves  upward  from  ground  to  cloud.  This  visual  display 

is  popularly  termed  "lightning."  The  return  stroke  actually  moves  downward, 

1 5 

typically  carrying  negative  charges  from  cloud  to  ground. 


11.  Horner,  F.  (1905)  Planet  Space  Sci.  13;1 137-1 150. 

12.  Huschke,  R.E.,  Ed.  ( 1 959)  Glossary  of  Meteorology,  Am.  Meteor.  Soc. , 

Boston,  038  pp. 

13.  Ishikawa,  H.  (1900)  Proc.  Res.  Inst.  Atmos.,  Nagoya  Enlv. , 8A:  1-274. 

14.  Pierce,  E.T.  et  al  (1902)  Final  Kept.  SRI  Project  No.  3738,  Contract  No. 

AF33(057)-7009,  Stanford  Research  Institute,  Menlo  Park,  Calif.,  132  pp. 

15.  Chalmers,  J.A.  (1907)  Atmospheric  Electricity,  2nd  ed. , Pergamon,  New 

York,  515  pp. 
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There  are  multiple  return  strokes  in  most  cloud -to-ground  discharges.  The 

, 1 <•  » on  16,  17,  18,  19 

number  of  return  strokes  may  vary  from  one  to  as  many  as  20. 

20 

Pierce  developed  empirical  rules  for  the  variation  with  latitude  of  the  average 
number  of  return  strokes  per  cloud-to-ground  discharge  and  the  proportion  of  all 
discharges  that  go  to  ground.  He  found  that  with  geographic  latitude  the  proportion 
of  discharges -to-ground  increases,  and  the  number  of  return  strokes  per  discharge 
decreases. 

21 

The  nature  of  the  cloud  discharge  is  not  clearly  understood.  Smith  observed 
that  most  (5(>  percent)  cloud  discharges  transfer  negative  charge  upward  from  a 

negative  center  of  charge  near  the  cloud  base  to  a positive  center  above.  Ogawa 
22 

and  Brook  observed  that  the  direction  of  propagation  of  charge  was  predominantly 

(7S  percent  of  cases  studied)  from  the  positive  center  of  charge  downward  to  the 

23 

negative  center  of  charge.  Takagi  observed  in  Japan  that  the  main  process  of 

cloud  discharge  involves  a positive  streamer  that  propagates  downward  from  a 

24 

positive  to  a negative  center  of  charge.  Moyer  observed  from  aloft  a line  storm 
in  which  about  75  percent  of  the  discharges  appeared  to  parallel  the  bases  of  single 
clouds  with  the  stroke  typically  emerging  from  one  extremity  and  re-entering  an 
opposite  extremity.  This  observation  was  made  on  12  November  1972  at  a location 
south  of  Dallas,  Texas. 

1.2.2  ATMOSPHERICS 

Atmospherics  are  electromagnetic  signals  which  emanate  from  lightning 

discharges.  The  term  "atmospherics"  is  frequently  shortened  to  "sferics." 

25 

Popov,  working  at  the  Pavlovsk  Magnetic  and  Meteorological  Observatory,  was 
the  first  investigator  to  study  the  phenomena  with  the  use  of  a detector.  Recent 
advances  in  sferics  work  have  been  largely  refinements  and  extensions  of  pioneer 

3 

work  accomplished  before  1940.  According  to  Horner,  the  principal  early 
workers  were  Appleton  and  Watson-Watt  in  the  United  Kingdom,  .Schonland  in 


16.  Pierce,  E.  T.  (1955)  Quart.  J.  Roy.  Meteorol.  .Soc.  ^:211-228. 

17.  MacKerras,  D.  (1968)  J.  Geophys.  Res.  ^:1175-1183. 

18.  Takeuti,  T.  (1965)  Proc.  Res.  Inst.  Atmos.,  Nagoya  University,  12A:l-70. 

19.  Workman,  E.  G. , Brook,  M, , and  Kitagawa,  N.  (1960)  J.  Geophys.  Res. 

65:1513-1517. 

20.  Pierce.  E.T,  (1970)  J.  Appl.  Meteor.  9:194-195. 

21.  .Smith,  L.G.  (1957)  Quart.  J.  Roy.  Meteor.  Soc.  83:103-111. 

22.  Ogawa,  T.,  and  Brook,  M.  (1964)  J.  Geophys.  Res.  69:5141-5150. 

23  Takagi,  M,  (1961)  Proc.  Res.  Inst.  Atmos..  Nagoya  University,  8B:1  -105. 

24.  Moyer,  V.  E.  (1974)  Personal  communication. 

25.  Popov,  A,.S.  (1896)  J.  Russ.  Phys.  Chem.  Soc.  28:7-9. 
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South  Africa,  Rureau  in  France,  Lugeon  in  Switzerland,  Austin  in  the  United 
States,  and  Norinder  in  Sweden. 

A common  method  used  to  study  sferlcs  is  the  examination  of  the  amplitude 
frequency  spectrum  of  the  received  sferlc  waveform.  A broadband  receiver  may 
be  used  to  record  incident  sferic  waveforms  for  analysis.  By  using  electromag- 
netic propagation  laws  for  very  low  frequency  (VI. F),  one  may  make  inferences 
concerning  the  amplitude  spectrum  of  the  source  from  a study  of  the  received 

. 14 

spectrum. 

I.eushin^*’  studied  the  geographic  correspondence  of  centers  of  sferic  activity 

and  thunderstorms  in  the  Soviet  Union  and  concluded  that  there  is  a one-to-one 

27 

correspondenci  between  them.  Barkalova  found  that  there  was  a good  correla- 
tion (ranging  from  0.  73  to  0.  97)  between  the  nu.nber  of  thunderstorms  and  the 
count  of  sferics  (recorded  above  a threshold)  at  two  stations  in  the  Soviet  Union. 

It  should  be  mentioned  that  sferlcs  may  originate  in  Sun  or  the  other  stars  as  well 
as  in  lightning. 

1.2.3  SFERICS  IN  THE  FORM  OF  DISCRETE  PUl.SES 

Sferics,  in  general,  are  continuous,  electromagnetic  phenomena.  Cloud-to- 
ground  discharges,  or  more  precisely,  the  associated  return  strokes,  produce 
intense  energy  in  the  VI. F portion  of  the  spectrum.  The  noise  thus  produced  is  of 

the  form  of  intermittent  pulses  superimposed  on  the  continuous,  background 
9ft2Q'^02'?  ‘11 

noise.  ’ ’ ’ Horner  and  Bradley'  report  that  the  transition  from  essen- 

tially discrete  pulses  to  continuous  waveforms  occurs  between  40  and  550  kHz. 

The  return  stroke  of  the  cloud-to-ground  discharge  also  generates  significant 
sferlcs  (K  pulses)  at  VLF  that  are  perhaps  one-tenth  the  magnitude  of  the  largest 
pulses.  Recoil  streamers,  which  occur  when  a propagating  streamer  meets  a 
center  of  charge  opposite  in  sign  to  the  streamer,  produce  relatively  small  K 
pulses.  The  leader  also  generates  K pulses,  but  these  are  very  small  in  magni- 
tude. Although  the  various  types  of  K pulses  are  much  less  intense  than  the 

large,  discrete,  VI. F pulses,  they  are  more  numerous  and  must  be  considered 

32 

when  analyzing  received  sferics, 

2fi.  I.eushin,  N.  E.  (1964)  AFCRI.,  Contract  AF19(628)-3880,  Am.  Meteor,  Soc. 
Translation. 

27.  Barkalova,  K.  N.  (1964)  AFCRL,  Contract  AF19(628)3880,  Am.  Meteor.  Soc. 

Translation. 

28.  Horner,  F.  (1958)  J.  Atmos.  Terr.  Phys.  13:140-154. 

29.  Malan,  D.J,  (1958)  Recent  Advances  in  Atmospheric  Electricity,  Pergamon, 

New  York,  557-56Ti 

30.  Kitagawa,  N.,  and  Brook,  M.  (1960)  J.  Geophys.  Res.  65:1927-1931. 

31.  Horner,  F. , and  Bradley,  P.A.  (1964)  J.  Atmos.  Terr.  Phys.  26:1155-1166. 

32.  Pierce,  E.T.  (1973)  Personal  communication. 


The  cloud  discharge  generates  quasl-diacrete  K pulses  at  VI. F that  are 

remarkably  similar  to  those  of  the  cloud-to-ground  discharges.  The  recoil 

14 

streamer  process  is  thought  to  cause  these  pulses.  It  is  clear  that  the  dominant 

29  33 

energy  at  VLF  generated  by  cloud  discharges  is  due  to  K pulses. 

In  summary,  both  cloud-to-ground  and  cloud  discharges  generate  significant 
sferics  at  VLF  in  the  form  of  quasi-discrete  pulses.  The  largest  pulses,  associ- 
ated with  cloud -to -ground  discharges,  are  due  to  the  initial  and  subsequent  return 
strokes.  Both  the  return  stroke  (cloud-to-ground  discharge)  and  recoil  streamer 
(cloud -to-ground  and  cloud  discharges)  produce  quasi-discrete  K pulses,  wfhlch 
are  less  intense  but  more  numerous  than  the  "main,"  VLF  pulses.  The  leader 
stage  of  a discharge  produces  sferic  pulses  at  VLF,  but  these  are  very  weak  in 
magnitude. 

1.3  llhjeclive  and  Scope  of  Report 

The  objective  of  this  study  was  to  contribute  information  on  the  following 
aspects  of  the  climatology  of  thunderstorms  and  lightning,  with  special  emphasis 
on  their  spatial  distribution:  an  estimation  of  the  distribution  over  much  of  the 
Eastern  Hemisphere,  thunderstorm  days;  the  number  of  lightning  discharges;  the 
incidence  of  lightning  discharge;  and  the  areal  concentration  of  lightning  discharge. 


i.  MIL  SKF.RH.S  l»\T  \ 

2.1  Dcfinilinn  and  Source 

The  sferics  data  which  were  available  for  this  study  are  counts  of  electromag- 
netic signals  (above  a threshold  intensity)  which  emanate  from  lightning  discharges. 
Personnel  of  the  U.  .S.  Air  Force  collected  these  data  in  1972  from  a large  network 
of  sferic  sensors  capable  of  recording  incident  sferics  on  a global  basis. 

2.2  The  Sferic!,  RccordcrK 

Although  sophisticated  electronically,  the  recorders  are  essentially  cathode- 
ray  tube  direction  finders  of  the  type  introduced  by  Watson-Watt  in  the  United 
Kingdom  more  than  40  years  ago.  The  recorders,  called  signal  monitors,  are 
passive  remote  sensors  that  detect,  classify,  and  record  the  vertically  polarized 
component  of  incident  sferics.  The  location  on  Earth  of  a sferic  is  determined  by 


*The  Eastern  Hemisphere  was  selected  because  of  the  availability  of  data  from 
that  area. 

33.  Horner,  F.  (1990)  J.  Atmos.  Terres.  Phys.  21:13-25. 


the  use  of  three  or  more  stations.  The  time  of  arrival  at  each  station  is  recorded: 
then  the  difference  in  time  of  arrival  between  stations  is  converted  to  distance  by 
the  use  of  the  estimated  propagation  rate  at  VLF.  Finally,  the  location  or  "fix"  of 
the  sferlc  is  determined  by  a straightforward  application  of  spherical  trigonometry. 
Sferics  incident  upon  the  antennas  are  recorded  on  a selective  basis.  The  criteria 
for  recording  a signal  are  the  azimuth  of  arrival  (sectoring),  the  waveform  cycle 
characteristics,  and  the  relative  intensity  of  the  waveform  (thresholding).  All 
incident  sferics  that  are  recorded  have  peak  energies  in  the  broadband  frequency 
range  of  250  Hz  to  60  kHz.  Most  of  the  energy  from  VLF  sferics  is  concentrated 
near  6 kHz,  ^ 

All  sferic  waveforms  recorded  by  the  signal  monitors  satisfied  the  time  and 
half-cycle  criteria  illustrated  in  Figure  1,  Also  in  Figure  1 are  examples  of 
actual,  large-source-strength  sferic  waveforms  taken  from  a paper  by  Kalakowsky 
and  Lewis.  The  waveforms  of  VLF  sferics  (at  distances  greater  than  about 
60  km)  are  of  the  form  of  sinusoidal  waves  that  build  to  a peak  and  then  decay  to 
zero.  The  criteria  for  recording  a waveform  are  that  the  half-cycle  of  maximum 
amplitude  must  occur  within  the  first  300  psec  of  the  waveform  and  that  the  total 
waveform  trace  must  be  no  greater  than  1000  psec  in  duration.  The  actual  wave- 
forms shown  in  Figure  1 were  recorded  simultaneously  at  Chicopee  and  Bedford, 
Massachusetts,  which  are  68.  6 mi  apart.  These  sferics  occurred  as  single  events 
associated,  usually,  with  cold  fronts  and  squall  lines,  in  close  proximity  to  the 
monitored  area.  Since  these  sferics  were  very  large  in  peak  field  strength 
(20  V m"^  at  50  nmi  minimum),  it  is  not  surprising  that  all  of  them  meet  the  cri- 
teria shown  in  Figure  1. 

2.3  \ri-a  of  Dais  Cvvrriigi* 

The  area  of  data  coverage  for  this  study  is  shown  in  Figure  2.  The  area  of 
coverage  was  determined  by  a combination  of  simple  azimuthal  sectoring  and  sys- 
tem thresholding.  The  best  estimate  of  the  area  of  data  coverage  in  the  mean  for 
1972  in  the  Eastern  Hemisphere  is  that  area  above  the  "threshold-sector  lines"  in 
Figure  2.  Data  generally  were  available  north  of  the  broken  line  in  Figure  2 for 
the  period  January  to  July  and  north  of  the  dash-dot  line  for  the  period  August  to 
December,  Although  the  available  data  covered  much  of  the  world,  exclusive  of 
parts  of  the  Americas,  only  the  Eastern  Hemisphere  was  considered. 

The  threshold  intensity,  in  decibels,  at  each  grid  point  is  a mean  estimate  of 
the  minimum  field  strength  below  which  no  signal  may  be  recorded.  Virtually  all 


34.  Bailey,  T.W.  (1972)  Personal  communication. 

35.  Kalakowsky,  C.B.,  and  Lewis,  E.A.  (1966)  Physical  Science  Research 

Papers.  No.  261.  AFCRL  Project  4603. 
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Figure  1.  The  Time  and  Half-cycle  Criteria  (a)  for  a Sferic 
Waveform  to  be  Recorded;  and  Examples  (b)  of  Intense  Wave- 
forms Recorded  at  (1)  Chicopee  and  (2)  Bedford,  Massachu- 
setts, from  the  Same  Source  Thunderstorms.  Arrows  in  (b) 
are  300  psec  apart  and  the  vertical  scale  is  in  V m"^.  (After 
Kalakowsky  and  Lewis,  1966) 

sferics  sensors  operate  above  a threshold  limit,  including  the  original  instrument 
designed  and  used  by  Popov  in  1896.  The  threshold  values  used  in  this  study  were 
determined  by  personnel  of  the  U.S.  Air  Force  using  a computer-based  model. 
The  model  inputs  were  the  geographic  location  of  the  sensors,  propagation  laws  at 
VLF,  and  the  individual  sensor  thresholds  (20  dB  minimum  above  a 1 mV  m * 
reference).  The  model  output  was  two  system  threshold  arrays  (over  the  grid 
shown  in  Figure  2)  that  give  mean  threshold  values  fur  the  periods  January  to  July 
and  August  to  December  of  1972.  Sample  calculations  using  these  arrays  will  be 
presented  in  Section  3. 
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Figure  2.  The  Area  of  Data  Coverage  in  the  Eastern  Hemisphere  with  Grid  Points. 
Each  grid  point  Is  centered  within  a data  block  bounded  by  parallels  and  meridians 
10-deg  apart.  The  broken  line  is  an  estimate  of  the  threshold-sector  line  for  the 
period  January  to  July  1972.  The  dash-dot  line  is  the  threshold-sector  line  for  the 
period  August  to  December  1972.  No  data  were  expected  south  of  the  threshold- 
sector  lines  during  the  respective  periods. 


2. 1 \niir»rKHli<>n  und  \olum<‘  of  Hula 

There  are  over  150,  000  datum  points  (individual  counts  of  signals)  available 
over  the  Eastern  Hemisphere  for  1972.  The  data  are  aggregated  in  four  synoptic 
time  steps  dally  (0000-0600Z.  0000-1200Z,  1200-1800Z,  1800-2400Z)  and  the  areal 
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resolution  U by  10-deg  latitude  and  longitude  blocks  centered  about  the  grid  points 
shown  in  Figure  2. 

A summary  of  the  completeness  of  the  available  data,  by  synoptic  time  incre- 
ments, is  given  in  Table  1.  Missing  data  occurred  either  at  random  or  toward  the 
middle  of  the  month.  Therefore,  straightforward,  linear  corrections  for  missing 
data  were  used  where  necessary. 


Table  1.  Summary  of  Number  of  Calendar  Days  Data  were  Available  by 
Synoptic  Time  Increments  for  Each  Month  of  1972  and  for  the  Year 


Month 

Days 

0000- 

0600Z 

0600- 

1200Z 

1200- 

1800Z 

1800- 

2400Z 

Jan. 

31 

28 

28 

27 

28 

Feb. 

29 

29 

29 

29 

29 

Mar. 

31 

31 

31 

31 

31 

Apr. 

30 

30 

30 

30 

30 

May 

31 

31 

31 

31 

31 

June 

30 

27 

27 

26 

26 

July 

31 

23 

24 

21 

23 

Aug. 

31 

31 

31 

31 

31 

Sep. 

30 

23 

23 

23 

23 

Oct. 

31 

29 

30 

28 

29 

Nov. 

30 

30 

30 

30 

30 

Dec. 

31 

30 

30 

31 

29 

Annual 

366 

342 

344 

338 

340 
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3. 1 Thunilef.lonii  Days 

An  estimate  of  the  number  of  thunderstorm  days  at  each  grid  p)oint  (with  avail- 
able data)  was  determined  for  selected  months  by  simply  counting  the  number  of 
days  on  which  signals  occurred.  The  IBM  360/65  computer  of  the  Data  Processing 
Center  was  programmed  to  count  the  days  with  signals.  A practical  constraint  was 
imposed:  that  the  sum  of  the  signals  for  one  day  must  be  at  least  three  in  the 
Northern  Hemisphere  to  be  counted.  Thus,  a zero  was  entered  at  those  grid  points 
where  the  sum  of  the  signal  count  for  one  day  (for  example,  the  sum  of  the  data  for 
0000-0600Z,  0600-1200Z,  1200-1800Z,  1800-2400Z)  was  less  than  three.  No 
similar  constraint  was  Imposed  for  the  Southern  Hemisphere  data. 
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The  reason  for  the  constraint  in  the  Northern  Hemisphere  was  that  north  of 
the  equator  the  threshold  values  were  sufficiently  low  allowing  several  spurious 
and  erroneous  signals  to  be  recorded.  An  important  feature  of  the  system  thresh- 
old arrays  was  that  relatively  few  signals  were  "masked  out"  in  north  latitudes, 
but  many  signals  were  masked  in  south  latitudes.  In  fact,  any  signal  recorded  in 

high  south  latitudes  (say,  poleward  of  45S)  was  of  relatively  great  field  strength 

34 

and  certainly  of  thunderstorm  origin.  ‘ The  effect  of  thresholding  will  be  illus- 
trated further  in  sample  calculations  in  this  section. 

The  mean,  global,  thunderstorm -day  chart  of  the  WMO  for  January,  is  pre- 
sented in  Figure  3.  The  analyzed  results  of  this  study  for  January  are  presented 
in  P'igure  4.  It  must  be  recalled  that  the  results  of  this  study  are  for  the  Eastern 
Hemisphere  only.  Similar  charts  are  presented  for  the  other  classic,  climatologic 
months  of  April,  July,  and  October.  In  each  of  the  charts  produced  in  this  study, 
there  appears  a heavy,  dash-dot  line  (threshold-sector  line)  which  delineates  the 
line  south  of  which  the  combination  of  azimuthal  sectoring  and  field-intensity 
thresholding  should  have  caused  a data  void.  Again,  recall  that  the  threshold  esti- 
mates were  averaged  over  two  periods:  January  to  July  and  August  to  December, 
1972.  Therefore,  the  extent  to  which  the  January  results  in  Figure  4 were  avail- 
able "south  of  the  line"  may  have  been  a qualitative  estimate  of  the  closeness  of 
the  actual  January  data  coverage  to  the  estimate  in  the  mean  for  the  period  January 
to  July.  On  the  other  hand,  the  amount  of  data  below  the  line  may  give  an  indica- 
tion of  the  sensitivity  of  thresholding.  That  is,  if  one  assumes  that  the  average 
threshold  array  for  the  period  January  to  July  was  very  close  to  that  for  January 
and  that  the  estimate  of  azimuthal  sectoring  was  accurate,  the  degree  of  "spillover" 
south  of  the  line  gives  a feel  for  the  error  in  the  threshold  array  estimate. 

It  is  very  likely  that  the  spill-over  effect  is  due  to  a combination  of  the  causes 
outlined  in  the  foregoing  paragraph.  Regardless  of  the  cause  or  causes  of  spillover, 
it  was  clear  that  the  recorded  data  were  valid;  the  problem  arises  in  their  treat- 
ment. 

For  the  present  purpose  of  estimating  the  thunderstorm -day  distribution  for 
selected  months,  the  spill-over  problem  is  not  very  significant,  because  the  analy- 
sis involves  only  the  counting  of  signals.  However,  the  problem  of  spillover  be- 
comes important  later  in  this  section.  The  results  which  appear  below  the  thresh- 
old-sector line  are  counts  of  days  on  which  Intense  VLF  signals  due  to  lightning 
discharges  (namely,  thunderstorms)  were  recorded.  The  only  ambiguity  arises  at 
grid  points  below  the  threshold -sector  line  where  blank  spaces  occur.  Here  one 
does  not  know  whether  this  marks  a true  data  void  (that  is,  no  signal  occurred),  or 

^Reproduced  from  the  Handbook  of  Geophysics  for  Air  Force  Designers  (1960). 


F'igure  3.  Thunderstorm-day  Chart  of  the  WMO  for  January  (all  years 
of  record)  for  the  World.  Isolines  are  in  thunderstorm  days.  (After 
Handbook  of  Geophysics,  1960) 


the  efrects  of  thresholding  and  sectoring  in  combination,  eliminating  signals  pro- 
duced by  lightning  discharges. 

The  major  thunderstorm  areas  in  the  world,  as  determined  by  a review  of  the 
WMO  data.  ^ are  in  North  and  South  America,  in  West,  Central,  and  South  Africa, 
and  in  Southeast  Asia.  These  areas  are  hereafter  termed  thunderstorm  foci 
because  they  dominate  the  world  distribution  of  thunderstorm  occurrence  on  a 
monthly,  seasonal,  and  yearly  basis. 

11.2  b\uluation  of  Krsbits 

When  comparing  the  WMO  thunderstorm-day  charts  to  those  of  this  study,  one 
should  keep  in  mind  that  both  have  inherent  advantages  and  disadvantages.  The 
primary  strength  of  the  former  is  that  they  cover  extensive  periods  of  record  at 
many  stations.  The  years  of  record  vary  widely,  however,  from  a low  of  1 yr  at 
Agedabia  in  Africa  (30N,  20E)  to  52  yr  at  Ufa  in  Russia  (54N,  45E).  The  main  weak- 
nesses of  the  WMO  effort,  as  cited  in  Section  1,  are  the  limitations  associated  with 
the  observational  technique  (the  hearing  and  recording  of  thunder)  and  the  extensive 
areas  over  the  globe  that  are  void  of  data.  The  very  low  c Dcurrence  of  thunder- 
storm activity  over  the  oceans,  reflected  in  the  WMO  compilation,  probably  is  due 
to  the  lack  of  observations. 

The  data  of  this  study  are  internally  complete  in  the  spatial  dimension,  but 
cover  only  1 yr  of  record.  The  vagaries  intrinsic  in  the  sensing  of  transient 
sferics  may  lead  one  to  question  the  efficacy  of  the  observational  techniques  used 
here.  However,  it  appears  that  the  remote  sensing  of  transient  sferics  is  a more 
sound  observational  approach  than  the  human  observation  of  thunder. 

A comparison  of  results  with  those  of  the  WMO  for  the  month  of  January  can 
be  made  by  looking  at  Figures  3 and  4.  Some  salient  observations,  by  geographic 
areas,  are: 

3.2.1  SOUTH  AFRICA 

The  correspondence  is  close,  in  order  of  magnitude,  between  the  South  African 
focus  of  the  WMO  and  that  of  this  study.  The  WMO  central  value  of  25  thunder- 
storm days  over  Madagascar  (15S,  45E)  in  Figure  3,  is  in  close  agreement  with  the 
value  of  28  in  Figure  4.  The  WMO  Isolines  terminate  sharply  east  of  Madagascar 
in  the  Indian  Ocean,  whereas  the  smoothed  analysis  of  this  study  indicates  a grad- 
ual decrease  of  activity  eastward,  stretching  across  the  ocean  to  India,  In  both 
studies,  the  activity  decreases  sharply  at  the  southern  tip  of  Africa  (ca.  35S). 

3.2.2  WE.ST  AFRICA 

Again,  there  is  agreement  in  the  order  of  magnitude  at  the  West  African  focus 
of  thunderstorm  days  between  the  two  studies.  The  value  of  this  study  in  Figure  4 
at  5N,  5E  of  ten  thunderstorm  days  is  close  to  the  Isoline  of  ten  in  the  WMO  chart. 
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Both  results  show  a drop  in  activity  to  the  immediate  north.  The  WMO  chart 
shows  a central  value  of  15,  whereas  the  maximum  value  of  this  study  is  ten  in 
West  Africa. 

3.2.5  NORTH  AFRICA,  THE  MKDITKRRANEAN  SEA,  AND  ARABIA 

There  are  no  thunderstorm  days  shown  for  January  over  North  Africa  in  the 
WMO  compilation  in  Figure  3.  A range  from  zero  to  11  thunderstorm  days  is 
indicated  by  this  analysis  in  Figure  4 for  the  same  area  (that  is,  roughly  from  15N 
to  30N  latitude  in  Africa).  The  central  value  in  the  Mediterranean  of  five  in  the 
WMO  chart  is  also  much  less  than  the  very  high  value  produced  in  this  study  of 
26  days  (35N,  15E)  with  at  least  one  thunderstorm.  The  most  remarkable  differ- 
ence that  shows  up  consistently  is  the  estimated  number  of  thunderstorm  days  over 
Arabia.  The  WMO  shows  no  thunderstorm  days  for  January  in  Figure  3,  whereas 
this  study  shows  a central  value  of  It)  at  25N,  45E  in  Figure  4. 

3.2.4  CONTINENTAL  El'RASIA 

In  general,  the  estimates  of  thunderstorm  days  produced  herein  are  higher 
across  Eurasia  than  those  of  the  WMO.  Both  results  show  a near  paucity  of 
thunderstorms  in  high  latitudes  of  Eurasia  for  January.  By  following  the  isoline 
of  one  thunderstorm  day  in  the  Northern  Hemisphere  in  both  charts,  one  may 
observe  that,  especially  in  Eastern  Eurasia,  the  results  of  this  study  indicate 
thunderstorm  occurrence  much  farther  poleward  than  does  the  WMO. 

3.2.5  INDIA 

The  estimates  of  the  WMO  and  this  study  of  one  or  more  thunderstorm  days  in 
Central  and  Southern  India  are  in  close  agreement.  As  mentioned  earlier,  the 
results  of  this  study  indicate  high  thunderstorm  activity  in  the  Indian  Ocean  south 
of  Ceylon,  whereas  the  compilation  in  the  mean  of  the  WMO  shows  no  days  with 
thunderstorms  south  of  5S.  The  WMO  value  of  one  thunderstorm  day  in  Northern 
India  (for  example,  at  Cherranpunji,  approximately  25N,  91E)  differs  from  the 
estimate  of  this  study  of  from  six  to  seven  thunderstorm  days  in  India  along  25N 
latitude  in  the  foothills  of  the  Himalayas. 

3.2.6  JAPAN 

A well-defined  center  of  thunderstorm  days  of  15  off  the  east  coast  of  Japan 
(35N,  145E)  is  estimated  herein  for  January,  whereas  the  WMO  chart,  shown  in 
Figure  3,  indicates  only  one  to  five  thunderstorm  days.  There  are  no  thunderstorm 
day.s  in  January  over  Kyushu,  Japan,  in  the  Yellow  Sea,  or  in  Korea  in  the  mean, 
according  to  the  WMO,  whereas  in  this  study  it  is  estimated  that  there  are  14  days 
at  35N,  135E  and  10  days  at  35N,  125E.  Moreover,  the  estimate  of  nine  thunder- 
storm days  In  the  North  Pacific  Ocean  along  35N  latitude  produced  in  this  study  is 
much  higher  than  the  estimate  of  one  thunderstorm  day  reported  by  the  WMO. 
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3.  2,  7 SOUTHEAST  ASIA 

There  id  no  generally  accepted  delineation  of  the  area  that  encompadses 
Southeadt  Asia.  For  purposes  of  discussion,  the  Southeast  Asian  focus  of  thunder- 
storm days  will  be  defined  arbitrarily  to  include  the  area  bounded  by  20N  and  lOS 
and  90E  and  175E.  The  estimate  in  this  study  of  one  thunderstorm  day  for  January 
in  the  Vietnam  area  (15N,  115E)  shown  in  Figure  4,  agrees  with  that  of  the  WMO 
shown  in  Figure  3.  However,  the  estimate  of  from  four  to  five  thunderstorm  days 
in  the  East  Indies  (for  example,  at  5S,  115E)  is  much  lower  than  the  WMO  estimate 
of  15  thunderstorm  days.  The  well-defined  isoline  of  five  thunderstorm  days  of 
this  study  shown  in  Figure  4 penetrates  southward  from  the  East  Indies  into  the 
Indian  and  even  Southern  Ocean  area.  The  isoline  of  five  thunderstorm  days  of  the 
WMO  does  not  penetrate  southward  of  about  15S  latitude  in  the  Indian  Ocean. 

3.2.8  AUSTRALIA 

There  is  close  agreement  in  the  vicinity  of  Australia  between  the  results  for 
January  of  the  W'MO  effort  and  this  study.  The  WMO  value  of  five  thunderstorm 
days  across  northern  Australia  (Figure  3)  compares  well  with  the  values  given  in 
this  study  (Figure  4),  ranging  from  two  to  five  thunderstorm  days  (that  is,  along 
15S  latitude  and  bounded  by  115E  and  155E  longitude).  The  value  of  roughly  five 
thunderstorm  days  in  the  area  of  the  Ellice  Islands  (ca  5S,  175E)  also  agrees 
closely  with  that  of  the  WMO. 

3.  2.  9 ANTARCTICA  AND  THE  SOUTHERN  OCEAN 

The  compilation  of  the  WMO  for  January  shown  in  Figure  3 indicates  no 
thunderstorm  days  south  of  50S.  In  this  study,  VLF  signals  of  the  type  which 
emanate  from  lightning  discharges  were  recorded  on  at  least  one  day  in  January 
throughout  much  of  the  Indian  and  Southern  Oceans  and  even  in  Antarctica,  as 
shown  in  Figure  4. 

The  WMO  compilations  of  thunderstorm  days  for  April  (Figure  5),  July 
(Figure  7),  and  October  (Figure  9)  were  compared  to  the  respective  estimates  of 
this  study  (P'igures  6,  8,  and  10)  in  a manner  analogous  to  that  given  above  for 
January.  A comparison  between  charts  by  classes  of  "five  thunderstorm  days" 

(that  is,  0-5,  5-10,  10-15,  and  on  on)  gave  a criterion  for  correspondence.  The 
correspondence  was  considered  "good"  if  values  were  within  a class  in  both  charts, 
"fair"  if  no  more  than  one  class  apart,  and  "poor"  if  more  than  one  class  apart. 

A summary  of  these  comparisons,  including  that  for  January,  is  given  in  Tables  2 
through  5. 

The  isoline  patterns  of  the  WMO  charts  and  those  of  this  study  differ  markedly. 
This  is  not  surprising  in  view  of  the  coarse,  grid  mesh  of  this  study  (that  is,  values 
by  10-deg  blocks  of  latitude  and  longitude)  and  the  concomitant  smoothed  analyses. 


< 

! 

i 

i 

I 


24 


P'igure  5.  Thunderstorm-day  Chart  of  the  V\'MO  for  April  (all  years  of 
record)  for  the  World.  Isolines  are  in  thunderstorm  days.  (After 
Handbook  of  Geophysics,  19fi0) 
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Figure  7.  Thunder.'jtorm-day  Chart  of  the  WMO  for  July  (all  years  of 
record)  for  the  World.  Isolines  are  in  thunderstorm  days.  (After 
Handbook  of  Geophysics,  1960) 
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Figure  8.  Thunderstorm-day  Chart  of  this  Study  for  July  1972  for  the 
Eastern  Hemisphere.  The  dash-dot  line  is  the  average,  threshold-sector 
line  for  the  period  January  to  July  1972.  Full  Isolines  are  in  thunder- 
storm days.  Broken  lines  are  intermediate  isoltnes 
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Fiffure  10.  Thunderstorm-day  Chart  of  this  Study  for  October  1972  for 
the  Eastern  Hemisphere.  The  dash-dot  line  is  the  average,  threshold- 
sector  line  for  the  period  August  to  December  1972.  Full  isolines  are 
in  thunderstorm  days.  Broken  lines  are  intermediate  isolines 
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Table  2.  Comparison  of  the  Compilation  of  Mean  Thunderstorm  Days  by  the  WMO 
over  the  Eastern  Hemisphere  for  January  and  the  Estimate  of  this  Study* 


Area 

Typical  Value  (or  Range) 
WMO  This  Study 

Correspondence 

1. 

South  Africa 
(south  of  equator) 

25 

28 

G 

2. 

West  and  Central 
Africa 

15 

10 

F 

3. 

North  Africa 
(north  of  1 5N) 

0 

O-ll 

P 

4. 

Mediterranean  Sea 

5 

20 

P 

5. 

A rabia 

0 

18 

P 

6. 

North  Eurasia 
(north  of  35N) 

0 

0-1 

G 

7. 

South  Eurasia 

1 

6 

F 

8. 

North  India 
(north  of  1 5N) 

1 

7 

F 

9. 

South  India 

1 

1 

G 

10. 

Japan 

1 

14 

P 

1 1. 

Southeast  Asia 

15 

5 

P 

12. 

North  Australia 
(north  of  25.S) 

5 

4 

G 

13. 

.South  Australia 

I 

1 

G 

14. 

Indian  Ocean 
(north  of  45S) 

1 

1-17 

P 

15. 

■Southern  Ocean 

0 

1-6 

P 

*A  judgment  of  the  correspondence  In  magnitude  of  the  estimates  is  given  as 
Good  (G),  Fair  (F),  or  Poor  (P)  as  explained  in  the  text.  The  Isoline  patterns 
differ  markedly. 
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Table  3.  Comparison  of  the  Compilation  of  Mean  Thunderstorm  Days  by  the  WMO 
over  the  Eastern  Hemisphere  for  April  and  the  Estimate  of  this  Study* 


A rea 

Typical  Value  (or  Range) 
WMO  This  Study 

Correspondence 

1.  South  Africa 

20 

21 

G 

(south  of  equator) 

2.  West  and  Central 

20 

7 

P 

Africa 

3.  North  Africa 

0,  >2 

1-5 

G 

(north  of  15N) 

4.  Mediterranean  Sea 

>2 

5 

G 

5.  Arabia 

0 

27 

P 

6.  North  Eurasia 

1 

1-11 

P 

(north  of  35N) 

7.  South  Eurasia 

1-10 

20 

P 

8.  North  India 

5 

28 

P 

(north  of  15N) 

9.  South  India 

5 

27 

P 

10.  Japan 

1 

7 

F 

11.  Southeast  Asia 

5-15 

1-29 

F 

12,  North  Australia 

1-5 

2-24 

P 

(north  of  25S) 

13.  South  Australia 

1-5 

2-22 

P 

14.  Indian  Ocean 

1-5 

9-24 

P 

(north  of  45S) 

15.  Southern  Ocean 

0-1 

1-15 

P 

♦A  judgment  of  the  correspondence  in  magnitude  of  the  estimates  is  given  as 
Good  (G),  Fair  (F),  or  Poor  (P)  as  explained  in  the  text.  The  isollne  patterns 
differ  markedly. 
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Table  4.  Comparison  of  the  Compilation  of  Mean  Thunderstorm  Days  by  the  WMO 
over  the  Eastern  Hemisphere  for  July  and  the  Estimate  of  this  Study 


Typical  Value  (or  Range) 

Area 

WMO 

This  .Study 

Correspondence 

1. 

South  Africa 
(south  of  equator) 

1-5 

1-5 

G 

2. 

West  and  Central 
Africa 

10-25 

2-16 

F 

3. 

North  Africa 
(north  of  15N) 

0-5 

0-3 

G 

4. 

Mediterranean  Sea 

0,  >2 

18 

P 

5. 

Arabia 

0 

0-7 

F 

6, 

North  Eurasia 
(north  of  35N) 

1-10 

1-23 

P 

7. 

South  Eurasia 

5-10 

23 

P 

8. 

North  India 
(north  of  15N) 

5-10 

13 

F 

9. 

South  India 

0-5 

6 

G 

10. 

Japan 

5 

14 

P 

11. 

.Southeast  Asia 

5-10 

0-15 

F 

12. 

North  Australia 
(north  of  25.S) 

1-5 

0-5 

G 

13. 

.South  Australia 

0,  >3 

0-4 

G 

14. 

Indian  Ocean 
(north  of  45.S) 

0,  >2 

0-3 

G 

15. 

■Southern  Ocean 

0-1 

0-3 

G 

♦A  judgment  of  the  correspondence  in  magnitude  of  the  estimates  is  given  as 
Good  (G),  Fair  (F),  or  Poor  (P)  as  explained  in  the  text.  The  isollne  patterns 
differ  markedly. 
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Table  5.  Comparison  of  the  Compilation  of  Mean  Thunderstorm  Days  by  the  WMO 
over  the  Eastern  Hemisphere  for  October  and  the  Estimate  of  this  Study 


Typical  Value  (or  Range) 

Area 

WMO 

This  Study 

Correspondence 

1. 

South  Africa 
(south  of  equator) 

5-15 

6-29 

P 

2. 

West  and  Central 
Africa 

5-25 

8-30 

G 

3. 

North  Africa 
(north  of  15N) 

0-5 

8-24 

P 

4. 

Mediterranean  Sea 

5 

21 

P 

5. 

Arabia 

0 

14 

P 

ti. 

North  Eurasia 
(north  of  35N) 

0-1 

1-21 

P 

7. 

South  Eurasia 

0,  >1 

12 

P 

8. 

North  India 
(north  of  15N) 

5 

11 

F 

9. 

South  India 

5 

30 

P 

10. 

Japan 

1 

12 

P 

11. 

Southeast  Asia 

5-10 

1-30 

P 

A judgment  of  the  correspondence  in  magnitude  of  the  estimates  is  given  as  Good 
(G),  Fair  (F),  or  Poor  (P)  as  explained  in  the  text.  The  isollne  patterns  differ 
markedly. 


together  with  the  extensive  areas  of  data  void  in  the  WMO  compilation.  There  also 
were  wide  differences  in  the  point  estimates  of  thunderstorm  days.  There  is  gen- 
erally good  agreement  in  the  comparison  of  the  results  for  the  thunderstorm-day 
focus  of  South  Africa  and  extremely  poor  agreement  In  Arabia.  There  were  16 
good,  10  fair,  and  30  poor  correspondences  of  estimates  shown  in  Tables  2 
through  5. 

One  could  ascribe  various  reasons  for  the  differences  in  correspondence  in 
thunderstorm  days  between  the  two  studies.  An  obvious  one  is  that  the  spatial  and 
temporal  differences  lead  one  to  expect  large  differences  in  results.  However, 
marked  differences  in  thunderstorm  days  at  a station  may  occur  from  one  year  to 
the  next,  due  to  many  factors,  including  the  vagaries  of  the  monsoon  regime. 
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Ramage^®  delineates  the  monsoon  area  as  that  between  35N  and  25.S  and  between 
30W  and  170E.  He  refers  to  a "bewilderingly  complex”  pattern  of  vertical  motion 
throughout  the  area  and  states  that  thunderstorms  typically  occur  rather  sporad- 
ically during  breaks  in  the  steady  regime  of  rain.  Outside  of  the  area  influenced 
by  the  monsoons,  there  also  were  considerable  differences  in  correspondence  of 
thunderstorm  days  between  the  two  studies.  It  seems  likely  that  many  of  the  differ- 
ences between  the  two  studies  as  outlined  herein  are  to  be  attributed  to  differences 
in  the  respective  space  and  time  considerations  of  the  observational  techniques. 

If  one  assumes  that  1972  was  not  an  anomalous  year,  then  it  is  clear  that  the 
great  difference  in  the  estimate  of  thunderstorm  days  for  Arabia  is  due  either  to  the 
lack  of  meteorological  stations  in  Arabia  (leading  to  a data  void  in  the  WIVIO 
data)  or  to  the  conclusion  that  the  sferics  recording  system  produced  erroneous 
counts  of  signals.  There  was  only  one  station  (at  Hail,  27N,  42E)  interior  to  Saudia 

Arabia  as  of  1953  (WMO,  1953).  The  other  11  stations  in  Saudia  Arabia  were  either 

37 

on  the  coast  or  on  adjacent  islands.  Griffiths,  based  on  years  of  experience  in 
Arabia,  stated  that  thunderstorms  may  occur  in  the  Yemeni  Highlands  in  Southwest 
Arabia  when  the  low-level  flow  is  northeasterly  from  the  Red  Sea  or  northwesterly 
from  the  Mediterranean  Sea. 

Both  thermally  induced  and  frontal  thunderstorms  are  quite  possible  in  desert 
areas,  although  the  probability  of  precipitation  is  slight.  An  analogy  is  that  it 
snows  quite  often  in  South  Texas,  but  snow  rarely  reaches  the  ground  as  precipita- 
tion. The  Pacific,  Indian,  and  Southern  Oceans  are  areas  where  the  differences 
between  the  two  studies  in  the  thunderstorm -day  distribution  were  almost  surely 

7 

due  to  a lack  of  data  in  the  WMO  compilation. 

3.3  The  Number  ni  l.ixhininx  Dixchurxex 

3.3.  1 THE  PHY.SICAL  MODEI. 

In  the  model  employed  here,  it  was  assumed  that  only  cloud  discharges  and 
cloud-to-ground  discharges  produce  significant,  intense,  VLF  sferics.  The  model 
is  illustrated  in  Figure  11. 

As  described  in  Section  1,  the  cloud  discharge  produces  quasi-discrete  K pulses 
in  the  VI. F range  of  the  spectrum.  The  cloud-to-ground  discharge  generates  the 
most  intense  (by  an  order  of  magnitude)  VI. F signals  as  well  as  weaker  K pulses. 
The  classic  cloud  model,  shown  in  Figure  11,  includes  a net,  negative  center  of 

charge  at  the  height  of  3 to  6 km  and  a net,  positive  center  of  charge  at  a height  of 
23  3 

7 to  11  km.  ■ ’ ■ Any  distribution  of  charge  would  be  acceptable  for  the  purpose  of 
this  study. 

36.  Ramage,  C..S.  (1971)  Monsoon  Meteorology,  Academic  Press,  New  York, 

296  pp. 

37.  Griffiths,  J.  F.  (1974)  Personal  communication. 
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Figure  11.  Model  of  the  Thundercloud  with  Electrical  Charge  Distribu- 
tion and  Types  of  Discharge.  The  latter  include  cloud  discharges  (1) 
and  (2)  and  cloud-to-ground  discharges  (3)  and  (4) 


3.3.2  THE  STATISTICAL  DISTRIBUTIONS 

The  first  step  in  the  estimation  of  the  distribution  of  lightning  discharges  was 
to  define  the  appropriate  statistical  distribution  of  the  received  quasi-discrete 
sferics.  Specifically,  the  distribution  of  received  sferics  due  to  K pulses,  the 
first  return  stroke,  and  subsequent  return  strokes  were  required.  The  selection 
of  the  appropriate  statistics  from  the  literature  is  presented  in  Table  6. 

There  is  wide  agreement  that  most  parameters  of  thunderstorm  activity  follow 
the  log-normal  probability  law.  Horner  and  Bradley  have  shown,  by  an  empirical 

O O 

study,  that  received  sferics  range  in  standard  deviation  from  6 to  8 dB.  Pierce 
states  that  K pulses  due  to  cloud  discharges  and  those  due  to  cloud -to -ground  dis- 
charges are  remarkably  similar  in  their  statistical  nature  due  to  the  observation 

38.  Pierce,  E,  T.  (1969)  Tech.  Rept.  2,  .SRI  Project  No.  7045,  Contract  No. 

DA.SA  01 -68-C-0073,  .Stanford  Res.  Inst.,  Menlo  Park,  Calif.,  90  pp. 
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Table  R.  Statistics  Used  to  Define  Distributions  of  .Sferlcs  Due  to  Return  Strokes 
and  K Pulses 


Statistics 

Value 

■Source 

Avg.  number  of  K pulses  per 
discharge  (cloud  or  cloud - 
to -ground) 

30 

(20-40  range) 

Pierce  (1969) 

Largest  field  strength 
intensity  of  K pulse 

0.  3 V m"^  at  100  km 
(49.  5 dB) 

Pierce  (1969) 

Mean  (Med.)  field  strength 
intensity  of  K pulse 

0.  06  V m ^ at  100  km 
(35.  5 dB) 

Estimated: 
ordinary  sta- 
tistical theory 

Mean  (Med. ) field  strength 
Intensity  of  first  return 
stroke  sferics 

3 V m at  100  km 
(69.  5 dB) 

Pierce  (1973) 
Horner (1964) 

Mean  (Med. ) field  strength 
intensity  of  subsequent 
return  stroke  sferics 

1.  5 V m ^ at  100  km 
(63.  5 dB) 

Pierce  (1969) 

Standard  deviation  of  field 
strength  intensity  of 
received  sferics  at  very 
low  frequency 

7 dB 

(6  to  8 dB) 

Horner  and 
Bradley  (1964) 

that  both  types  of  K pulse  are  caused  by  recoil  streamers,  as  explained  in  Section  1. 
However,  very  little  is  known  about  K pulses  in  relation  to  the  more  intense  types 
of  VLF  sferics. 

The  statistical  distributions  used  in  this  study  were  assumed  to  follow  the  log- 
normal law;  that  is,  the  peak  field  strength  of  all  received  sferlcs  follow  a normal 
probability  distribution  about  the  mean  (also  the  median)  with  values  expressed  in 
decibels  namely,  on  a logarithmic  scale.  The  distributions  are  defined  in  Table  7. 

The  values  of  field-strength  intensity  in  decibels  (dB)  given  in  Table  6 were 
determined  by  the  following  convention: 

dB  * 20  log  , (1) 

where  E is  the  voltage  of  the  peak  received  field  strength  of  the  waveform  above 
P .1 

the  reference  field  strength,  E^,,  of  1 mV  m . This  convention  was  also  used  in 
the  threshold  array  model.  The  mean  field  strength  value  for  the  K pulse  distribu- 
tion in  Table  6 was  calculated  by  assuming,  as  is  frequently  done,  that  the  largest 
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Table  7.  Definition  of  the  Statistical  Distributions  Used  in  this  Study,  (Values  are 
in  decibels.  Standard  deviation  is  denoted  by  SD) 


Sferic  Generator 

Distribution 

Mean 

•SD 

♦ 1 SD 

+ 2 SD 

+3  SD 

First  return 
stroke 

log-normal 

69.  5 

■ 

76.  5 

83.  5 

90.  5 

Subsequent  return 
strokes 

log-normal 

63.  5 

D 

70.  5 

77.  5 

84,  5 

K pulses 

log-normal 

35.  5 

D 

42.  5 

49.  5 

56,  6 

K pulse  intensity  (field  strength  of  0,3  V m”^)  occurs  two  standard  deviations 
from  the  mean.  Then,  with  a standard  deviation  of  7 dB,  one  obtains  a value  of 
35.  5 dB  (or  0.  06  V m"^)  for  the  mean  intensity  by  a straightforward  application 
of  the  normal  law. 

3.3.3  LATITUDINAL  VARIATION  OF  LIGHTNING  PARAMETERS 
20 

Pierce  gave  the  only  available  estimate,  based  on  empirical  studies,  of  the 
variation  with  latitude  of  the  average  number  of  return  strokes  per  cloud-to-ground 
discharge  and  the  proportion  of  all  discharges  that  go  to  ground.  The  line  of  best 
fit,  based  on  nine  sets  of  data,  for  the  average  number  of  return  strokes  per  cloud- 
to-ground  discharge,  n,  has  the  equation 

n * 6 - (ij./30)  , (2) 

where  ^ is  latitude  in  degrees.  The  equation  for  the  proportion  of  ground  dis- 
charges, P,  is 

P • 0.  1 ( 1 > ((i,/30)^l  , (3) 

where  ^ is  the  latitude  in  degrees.  This  equation  is  the  best  fit  for  19  sets  of  data. 
20 

As  Pierce  points  out,  other  factors  (orography,  type  of  thunderstorm,  and  so  on) 

ideally  should  be  considered  to  refine  the  "rough  fit"  represented  by  these  equations. 
3 8 

One  may  develop  from  Eq.  (3)  an  expression  for  the  ratio,  M,  of  the  number 
of  cloud  discharges,  D^,,  to  the  number  of  cloud-to-ground  discharges,  Dj^,  for  a 
given  latitude, 
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(4) 


M = 


Dc 


D, 


G 


Then 


, 9 - (^/30)f,  (5) 

1 + («/30)^ 

Since,  by  definition,  there  is  only  one  first  return  stroke  per  cloud -to -ground 
discharge,  it  follows  from  Eq.  (2)  that  the  number  of  subsequent  return  strokes 
per  discharge,  q,  for  a given  latitude,  <)>,  may  be  written 

q = 5 - (di/aO)  . <9) 

The  relationship,  in  the  mean,  among  the  first  return  stroke,  Rj,  the  subsequent 
return  strokes,  Rj,,  and  K pulses  (due  to  the  cloud  discharge  or  cloud -to-ground 
discharge)  is 

Rj:Rg:K  = l:q:30  , (7) 

where  q is  taken  from  Eq.  (6)  and  the  mean  value,  30,  from  Table  6. 

A sense  of  the  latitudinal  variation  of  the  lightning  parameters  given  by  Eqs. 

(2),  (3),  (5).  and  (6)  is  given  in  Table  8.  With  the  exception  of  the  proportion  of 

all  discharges  that  go  to  ground,  P,  which  increases  poleward,  all  parameters  in 

20 

Table  8 decrease  in  magnitude  from  equator  to  pole.  Pierce  stated  that  this 
increase  may  be  explained  by  a meteorological  argument:  The  smaller  the  dis- 
charge distance  in  the  cloud  (that  is,  from  the  net,  negative  center  of  charge  to  the 
net,  positive  center  of  charge  in  Figure  11)  relative  to  the  discharge  distance  from 
the  negative  cloud  center  of  charge  to  the  ground,  the  greater  the  probability  of  a 
cloud  discharge  relative  to  a discharge  to  ground.  This  concept  leads  one  to  infer 
that  the  height  of  the  net,  negative  center  of  charge  is  greater  on  the  average  in  the 
tropics  than  in  temperate  and  polar  latitudes,  although  this  has  not  been  demon- 
strated. It  should  be  mentioned  that  the  ratio  of  cloud  to  cloud -to -ground  dis- 
charges, M,  was  constrained  to  be  no  less  than  unity.  It  seems  most  unlikely  that 
there  are  more  ground  than  cloud  discharges,  on  the  average,  even  at  high 
latitudes. 
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Table  8.  Latitudinal  Variation  of  Lightning  Parameters  by  5-deg  Lati- 
tude Increments  from  Equator  to  Pole.  (The  values  are  applicable  in 
both  the  Northern  and  .Southern  Hemispheres) 


Latitude  (deg) 

n 

q 

P 

M 

0 

6.00 

5.  00 

0.  10 

9.  00 

5 

5.  75 

4.  75 

0.  10 

8.  73 

10 

5.  50 

4.  50 

0.  11 

8.00 

15 

5.  25 

4.25 

0.  12 

7.  00 

20 

5.  00 

4.00 

0.  14 

5.  92 

25 

4.  75 

3.  75 

0.  16 

4.  90 

30 

4.  50 

3.  50 

0.  19 

4.  00 

35 

4.  25 

3.  25 

0.  22 

3.  24 

40 

4.  00 

3.00 

0.  26 

2.  60 

45 

3.  75 

2.  75 

0.30 

2.08 

50 

3.  50 

2.  50 

0.35 

1.  65 

55 

3.  25 

2.  25 

0.  40 

1.  29 

fiO 

3.00 

2.00 

0.46 

1.00 

65 

2.  75 

1.  75 

0,  52 

1.  00 

70 

2.  50 

1.  50 

0.  59 

1.  00 

75 

2.  25 

1.25 

0.  66 

1.  00 

80 

2.00 

1.00 

0.  74 

1.00 

85 

1.  75 

0.  75 

0.  82 

1.00 

90 

1.  50 

0.  50 

0.  91 

1.00 

n Average  number  of  return  strokes  (first  and  subsequent)  per  cloud- 
to-ground  discharge. 

q Average  number  of  subsequent  return  strokes  per  cloud-to-ground 
discharge  (n-1). 

P Proportion  of  all  discharges  that  go  to  ground. 

M Ratio  of  the  number  of  cloud  discharges  to  the  number  of  cloud-to- 
ground  discharges. 
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3.3.4  DETEHMINATION  OF  THE  NFMBEK  OF  LIGHTNING  DISCHAHGE.S 

Two  procedures  for  estimating  the  number  of  lightning  discharges  at  each 
grid  point  were  developed  under  different  assumptions.  In  the  first  procedure. 
Case  1,  the  assumption  was  made  that  the  values  of  threshold  intensity  were 
sufficiently  large  to  preclude  the  sensing  of  K pulses;  that  is,  it  was  assumed  that 
only  sferics  due  to  the  intense  return  stroke  of  the  cloud-to-ground  discharge  were 
recorded  by  the  sferics  network.  The  absolute  minimum  threshold  value  in  the 
threshold  arrays  was  49.4  dB.  This  value  represents  that  intensity  in  field 
strength  below  which  no  signal  theoretically  could  have  been  recorded.  All  other 
threshold  values  were  50  dB  or  greater.  The  minimum  threshold  intensity  of 
49.4  dB  was  nearly  two  standard  deviations  greater  than  the  mean  of  the  K distri- 
bution, as  seen  from  Table  7.  Therefore,  if  the  sample  statistics  given  in  Table  G 
were  close  approximations  to  the  state-of-nature,  and  if  the  threshold  values  were 
reasonable  estimates  of  system  response  in  the  mean,  then,  clearly,  few  K pulses 
were  recorded  by  the  sferics  system. 

In  the  second  procedure.  Case  2,  the  assumption  was  made  that  K pulse 
sferics  were  recorded  as  well  as  sferics  due  to  return  strokes.  Case  1 was  used 
in  the  actual  analyses  of  the  data. 

In  Case  1,  we  express  the  number  of  VLF  signals  recorded,  X,  at  a given 
grid  point  as 

X = QjR,  r Q.,Rj.  , (8) 

where  Rj  is  the  first  return  stroke  (always  one  per  cloud-to-ground  discharge)  and 
R^  is  the  number  of  subsequent  return  strokes  per  discharge  (varies  with  latitude). 
The  symbols  Qj  and  Q,,  are  the  respective  probabilities  of  occurrence  of  a single 

R,  or  Rj,  sferic  above  a lower  limit  of  threshold  intensity. 

* ^ 39 

one  may  understand  better  Eq.  (8)  by  the  use  of  relative  probability  curves. 

Relative  probability,  as  used  here,  is  defined  as  the  probability  of  occurrence  of  a 
single  even  (Rj  or  R,,)  with  respect  to  a reference  (for  example,  the  mean,  field- 
strength  Intensity  of  the  Rj  distribution).  The  log-normal  curves  of  relative  prob- 
ability for  Rj  and  R^,  that  are  valid  along  the  equator  are  shown  in  Figure  12. 

(One  must  construct  a different  set  of  curves  for  each  of  the  18  latitude  points  in 
the  grid  in  Figure  2,  ) A given  threshold  intensity  serves  as  the  lower  limit,  in 
decibels,  of  the  probability  integral  for  a given  distribution.  For  example,  in  the 
typical  case  of  a threshold  intensity  of  .53.  G dB  at  55N,  35E,  the  probability  of 
occurrence  of  a first  return  stroke  greater  than  or  equal  to  53.  (>  dB  may  be 
expressed  according  to  the  log-normal  law  as 

39.  Pierce,  E.T.  (1968)  Tech.  Rept.  49,  SRI  Project  No.  4240,  Contract  No. 

DA  36-039-AMC:-00040E.  Stanford  Res.  Inst.,  Menlo  Park,  Calif.,  94  pp. 
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Intensity  of  signal  from  a cloud-to-ground  discharge 


Figure  12.  Probability  of  Occurrence  of  a Signal  above  a I hresh- 
old  Intensity  at  the  Ktjuator  Helatlve  to  the  Mean  (set  to  zero)  of 
First  Heturn  Stroke  Hist ribution  (solid  curve).  Broken  curve  is 
subsequent  return  stroke  distribution.  Latitude  is 


-T-/ 


Lp  - aL4L^ 


where  o i;i  th*-  standard  deviation  (assumed  to  be  7 d}3)  and  T,  in  decibels.  Is  the 
threshold  intensity.  The  curves  in  Figure  12  were  constructed  relative  to  the  mean 
of  the  first  return  stroke  distribution  (set  to  zero)  along  the  abscissa  and  in  accord- 
ance with  the  proportions  given  by  Kq.  (7)  for  the  ordinate  values.  Abscissa  values. 


therefore,  represent  decibel  values  "up"  or  "down"  the  logarithmic  scale  of 
Intensity  relative  to  the  mean  of  the  distribution  of  the  first  return  stroke. 

We  may  solve  for  Rj,  in  Eq.  (7)  in  terms  of  Rj  and  substitute  for  Rg  in  Eq.  (8) 
to  obtain 

X = Q,Rj  qQj,R,  . (10) 

Then  we  may  solve  Eq.  (10)  for  Rj,  which  is  equivalent  to  solving  for  the  number 
of  cloud-to-ground  discharges,  since  there  is  one  "first  return  stroke,"  Rj, 

per  D^: 


X 


(11) 


Now  by  the  substitution  of  the  appropriate  probability  integrals  in  the  denominator 
of  Eq.  (11),  Qj  and  Q^,,  with  the  concomitant  value  of  threshold  intensity  at  a grid 
point  as  the  lower  limit  in  each  integral,  we  may  solve  Eq.  (11)  for  at  the  grid 
point. 

The  final  step  in  Case  1 is  to  find  an  estimate  of  the  total  number  of  dis- 
charges, D,j,,  comprised  of  cloud  discharges,  D^.  and  cloud-to-ground  discharges, 
Dj^.  By  the  substitution  for  from  Eq.  (4)  we  have 

D_  D„  (1  * M)  . (12) 

i 

A sample  calculation  for  Case  1 of  the  estimated  number  of  lightning  dis- 
charges, D,p,  that  occurred  during  January  1972  within  the  data  block  centered  on 
25N,  125E  (1,  119,  165  km^)  is 

= 25N  (latitude) 

X = 688  (sum  of  all  VI. K signals  recorded  .at  25N,  125E  for  January  1972) 

X = '(69  (x  adjusted  for  missing  data  . . . the  average  data  coverage  for 
January  was  27.  75  from  Table  1 ...  a linear  adjustment  is 
(31/27.75)  (688)  or  768.6  signals) 

T - 66.9  dB  (threshold  value  at  25N,  125E  which  is  the  lower  limit  on  the 
probability  integrals) 

q « 3.75  (number  of  return  strokes  pei  cloud-to-ground  discharge  at  25N 
or  25S  latitude  . . . from  Table  8) 

M * 4,90  (ratio  of  the  number  of  cloud  to  cloud-to-ground  discharges  at  25N 
or  25S  latitude  . . . from  Table  8) 
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Then  the  procedure  used  to  obtain  Eq.  (9)  yields  a Qj  of  0.  944  and  a of  0,  316. 
(The  mean  of  the  first  return  stroke  distribution  is  69.  5 dll  from  Table  7.  There- 
fore, the  threshold  value  of  66.  9 dH  at  25N,  125E  is  ”2.  6 dB  down”  from  the  mean. 
One  may  estimate  Qj  and  by  looking  at  the  respective  areas  under  the  curves 
to  the  right  of  the  -2.  6 dB  abscissa  value  in  Figure  12).  We  now  may  solve  Eq. 

(12)  for  the  estimated  number  of  lightning  discharges,  U.j,,  that  occurred  during 
January  1972  within  the  data  block  centered  on  25N,  125E. 

769 

”t  o7644  * (3.  75)(0.  316)  . 

2480  . 

This  value  appears  in  Figure  13  at  25N,  125E,  rounded  off  to  the  nearest  100 
lightning  discharges  (that  is,  25). 

The  IBM  360/65  computer  was  programmed  to  solve  Eq.  (12)  for  the  estimated 
number  of  lightning  discharges,  D.^,,  at  each  grid  point  in  Figure  2 in  the  following 
expanded  form 


D 


T 


X 


dT  * q 


[1  + M|  , (13) 

where  the  unprimed  quantities  refer  to  the  distribution  of  the  first  return  stroke 
and  the  primed  quantities  refer  to  tha  of  the  subsequent  return  stroke;  X is  the 
number  of  VI, F signals  recorded;  a is  the  standard  deviation  (7  dB  in  both  distribu- 
tions); T^  is  the  threshold  value  at  a given  grid  point;  T is  the  intensity  of  the  mean 
(median)  field  strength  of  received  signals  in  decibels;  q is  the  number  of  subse- 
quent return  strokes  per  cloud-to-ground  discharge  given  by  Eq.  (6);  and  M is  the 
ratio  of  cloud  to  cloud-to-ground  discharges  given  by  Eq.  (5). 

In  Case  2,  we  assume  that  K pulses,  which  may  originate  in  both  cloud  dis- 
charges and  cloud-to-ground  discharges,  were  recorded  by  the  system.  Equation  (8) 
now  becomes 

X * Q,Hj  4 , (14) 

(1)  (2) 
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Term  (0  is  the  contribution  due  to  the  cloud-to-ground  discharge.  The  symbols 
Q,,  Q „ and  Q,.  are  the  probabilities  of  occurrence  of  the  first  return  stroke, 
subsequent  return  strokes,  and  K pulses,  respectively.  We  use  fij  and  as 
before,  and  K„  is  the  number  of  recorded  K pulses  due  to  the  ground  discharge. 

VJ 

Term  (2)  is  the  contribution  due  to  the  cloud  discharge.  Since  we  assumed  that 
the  statistics  of  K pulses  were  identical  regardless  of  origin,  is  the  same  as 
in  term  (1).  The  number  of  recorded  K pulses  due  to  the  cloud  discharge  is  repre- 
sented by 

Bv  substitution  for  R^,  K,,,  and  K-,  in  Eq.  (14)  with  their  equivalents  from 
Eq.  (7),  we  may  solve  for  the  number  of  cloud-to-ground  discharges,  D^.  (same  as 
H j ),  as  in  Case  1 : 


° Q,  + 3(Xij,.  + 30MQ|^  • 

A sample  calculation  for  Case  2,  which  is  hypothetical  since  Case  2 was  not 
used  in  the  analysis,  is  as  follows. 

I.et 

<l>  fiSN  (latitude) 

X = f)00  (sum  of  all  VI. F signals  recorded  during  one  month  at  a grid  point) 

T = 50.0  dB  (threshold  intensity) 

Following  the  procedure  in  Case  1,  we  obtain  these  values: 

Qj  = 0.99 

= 0.  97 

» 0.  02 

K, 

q = 1.  75  (from  Table  8) 

M 1.  00  (from  Table  8) 

Then  Eq.  (15)  may  be  solved  for  the  estimated  number  of  lightning  discharges,  D.^,; 

500 

•^T  “ 0.99  t (1.75)(0.97)  + 30(0.02)  + 30(1.0)(0.2)  • 

" 330  . 

The  result  for  D.j,  If  we  use  the  Initial  values  given  here,  but  under  the  Case  1 
assumption  that  no  K pulses  are  recorded,  is 
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= (186)(1  < 1.0)  = 372 


That  ia,  for  a grid  point  at  65N  with  an  associated  threshold  Intensity  of  50,0  dB, 
the  Case  1 procedure  leads  to  an  estimate  of  the  number  of  lightning  dlachargee, 
roughly  11  percent  higher  than  that  for  Case  2. 

•Since  the  very  low  threshold  intensity  of  50.  0 dB  leads  to  a difference  of  only 
11  percent  between  the  results  of  the  Case  1 and  Case  2 procedures,  it  was  evident 
that  the  choice  of  procedure  was  not  critical  to  the  work  as  long  as  the  threshold 
values  were  valid.  If  the  threshold  values  were  significantly  lower,  the  choice  of 
procedure  would  have  been  very  critical.  As  mentioned  earlier,  the  lowest  thresh- 
old intensity  at  a grid  point  was  49.4  dB.  All  other  values  were  greater  than 
50.0  dB. 

The  results  of  the  Case  1 type  analysis  of  the  sferlcs  data  are  given  in 

Figures  13  through  25  for  each  month  in  1972  and  for  the  year  as  a whole.  The 

background  map  shown  in  Figure  2 appears  in  each  of  these  figures.  .Smoothed 

2 6 

isolines  of  the  grid  point  results  are  given  in  powers  of  10  (that  is,  10  to  10  ), 
and  the  chart  entries  are  in  hundreds  of  discharges.  The  threshold-sector  lines 
that  appeared  in  the  thunderstorm-day  charts  (Figures  6 through  10)  were  repro- 
duced here  for  each  month.  It  is  to  be  recalled  that  this  line  for  a given  month  is 
an  estimate,  based  on  arrays  of  mean  threshold  intensity,  that  delineates  the  areal 
extent  of  data  coverage  (north  of  the  line)  due  to  the  combination  of  minimum  limits 
of  threshold  intensity  and  azimuthal  sectoring. 

Values  and  isolines  that  appear  south  of  the  threshold-sector  line  should  be 
viewed  with  extreme  caution.  The  extent  to  which  the  values  and  isolines  "spill 
over"  the  line  may  be  a measure  of  the  incorrectness  of  the  threshold  array  itself, 
or  it  may  simply  be  a result  of  individual  monthly  departures  from  the  mean  esti- 
mates. Due  to  the  practical  problem  that  no  values  of  threshold  intensity  were 
available  south  of  the  threshold-sector  line,  the  convention  was  adopted  that  the 
threshold  value  of  the  "three-standard-deviation  point"  of  each  of  the  statistical 
distributions  (given  in  Table  7)  was  used  everywhere  south  of  the  line.  The  pru- 
dent decision,  in  view  of  all  the  uncertainties,  may  be  to  ignore  all  results  below 
the  threshold-sector  line. 

It  is  further  recommended  that  all  values  that  appear  at  and  south  of  45S  lati- 
tude in  Figures  13  through  25  be  considered  questionable  until  confirmation  or 
rejection  of  them,  based  on  other  studies,  is  possible.  This  is  due  to  the  fact  that 
values  of  threshold  intensity  were  very  close  to  the  three-standard-deviation  point 
of  the  respective  distributions  given  in  Table  7 at  and  south  of  45S.  The  result  is 
that  only  a few  recorded  signals  south  of  45S  were  converted  to  many  discharges 
by  the  use  of  the  probability  integral  technique  explained  previously.  This  may 
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Figure  13.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for 
January  1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  pe- 
riod January  to  July  1972.  Values  are  In  hundreds  of  discharges  for  10-deg  data 
blocks  centered  at  the  points  shown.  Isolines  are  number  of  discharges  in  powers 
f 10 
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Figure  14.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for 
February  1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  pe- 
riod January  to  July  1972.  Values  are  in  hundreds  of  discharges  for  10-deg  data 
blocks  centered  at  the  points  shown.  Isolines  are  number  of  discharges  in  powers 
of  10 


Figure  15.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for 
March  1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  period 
January  to  July  1972.  Values  are  in  hundreds  of  discharges  for  10-deg  data  blocks 
centered  at  the  points  shown.  Isolines  are  number  of  discharges  in  powers  of  10 
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Figure  16.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for  April 
1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  period  January 
to  July  1972.  Values  are  in  hundreds  of  discharges  for  10-deg  data  blocks  centered 
at  the  points  shown.  Isolines  are  number  of  discharges  in  powers  of  10 
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Figure  17.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for  May 
1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  period  Janu- 
ary to  July  1972.  Values  are  in  hundreds  of  discharges  for  10-deg  data  blocks 
centered  at  the  points  shown.  Isollnes  are  number  of  discharges  in  powers  of  10 


Figure  18,  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for  June 
1972,  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  period  Janu- 
ary to  July  1972,  Values  are  in  hundreds  of  discharges  for  10-deg  data  blocks  cen- 
tered at  the  points  shown.  Isolines  are  number  of  discharges  in  powers  of  10 


Figure  19.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for  July 
1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  period  Janu- 
ary to  July  1972.  Values  are  in  hundreds  of  discharges  for  10-deg  data  blocks 
centered  at  the  points  shown.  Isollnes  are  number  of  discharges  in  powers  of  10 
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Figure  20.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for 
August  1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  period 
August  to  December  1972.  Values  are  in  hundreds  of  discharges  for  10-deg  data 
blocks  centered  at  the  points  shown.  Isolines  are  number  of  discharges  in  powers 
of  10 


Figure  21,  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for 
September  1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the 
period  August  to  December  1972.  Values  are  in  hundreds  of  discharges  for  10-deg 
data  blocks  centered  at  the  points  shown.  Isolines  are  number  of  discharj^es  in 
powers  of  10 


Figure  22.  Distribution  of  the  Estimated  Number  of  Lightning  Discharges  for 
October  1972.  The  dash-dot  line  is  the  average  threshold-sector  line  for  the  pe- 
riod August  to  December  1972.  Values  are  in  hundreds  of  discharges  for  10-deg 
data  blocks  centered  at  the  points  shown.  Isolines  are  number  of  discharges  in 
powers  of  10 
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Figure  23.  Diatributlon  of  the  Fatimated  Number  of  Lightning  Ulachargea  for 
November  1972.  The  ciaah-dot  line  ia  the  average  threahold -sector  line  for  the 
period  August  to  December  1972.  Values  are  in  hundreds  of  discharges  for  10- 
deg  data  blocks  centered  at  the  points  shown,  laolines  are  number  of  discharges 
in  powers  of  10 


F’ijjurp  25.  Oistribution  of  the  Estimated  Number  of  I.ightning  Discharges  for  1972 
The  broken  line  is  the  average  threshold-sector  line  for  the  period  January  to  July 
1972  and  the  daih-dot  line  is  that  for  the  period  August  to  December  1972.  Values 
are  in  hundreds  of  discharges  for  10-deg  data  blocks  centered  at  the  points  shown. 
Isolines  are  number  of  discharges  in  powers  of  10 


1 

I have  led  to  an  overestimation  of  the  number  of  lightning  discharges  if  the  received 

sferics  did  not  follow  the  log-normal  distribution  in  the  extreme  tails  of  the 
distribution. 

Even  a cursory  examination  of  the  threshold-sector  line  and  the  associated 
"spillover"  for  each  month  in  Figures  13  through  25  indicates  that  the  threshold- 
sector  estimate  was  quite  good,  indeed,  at  least  in  terms  of  its  effect  on  data 
coverage,  for  each  month  except  August  (Figure  20),  September  (Figure  21),  and 
I ■ October  (Figure  22).  In  these  months  there  was  considerable  spillover  of  data 

and  results  south  of  the  threshold-sector  line  and  east  of  95E  longitude.  It  seems 
likely  that  the  actual  threshold-sector  lines  of  the  system  for  these  months  were 
simply  departures  from  the  mean  estimate,  even  though  it  is  impossible  to  state 
this  with  any  assurance. 

It  should  be  pointed  out  that  the  combination  of  the  round-off  convention  used 
in  the  estimate  of  the  number  of  lightning  discharges  (that  is,  to  the  nearest  100 
discharges),  together  with  the  constraint  used  in  the  Northern  Hemisphere  for  the 
thunderstorm-day  estimates  (that  is,  at  least  3 VLF  signals  in  order  to  be  counted), 
resulted  in  two  monthly  charts  with  estimated  thunderstorm  days  but  no  estimated 
lightning  discharges  at  one  or  two  grid  points.  This  occurred  in  January  at 
35N,  105E  and  in  April  at  65N,  35E  and  at  45N,  45E, 

The  thunderstorm -day  charts  (Figures  6 through  10)  are  not  necessarily 
\ isomorphic  to  the  respective  lightning  discharge  charts  (Figures  13,  Ifi,  19,  and 

22)  in  every  respect.  The  thunderstorm  day  is  the  report  of  the  occurrence  of  at 
least  one  thunderstorm  during  a calendar  day  regardless  of  the  number  of  storms 
that  actually  occurred.  Therefore,  if  many  storms  occurred  on  a few  days  of  the 
month,  the  respective  distributions  of  thunderstorm  days  and  lightning  discharges 
could  be  quite  different. 

The  validity  of  the  estimates  of  the  number  of  lightning  discharges  given  in 
Figures  13  through  25  is  not  known,  although  the  values  of  the  incidence  of  dis- 
charge (to  be  given  later  in  Section  3)  which  are  based  on  the  estimates  of  the 
number  of  lightning  discharges,  appear  reasonable.  It  would  be  surprising  if  it 
were  demonstrated  that  the  distribution  of  the  number  of  lightning  discharges  had 
no  relationship  to  the  state-of-nature  even  though  the  magnitude  of  the  values  may 
be  doubtful;  that  is,  at  the  very  least,  it  seems  likely  that  the  results  presented  in 
Figures  13  through  25  represent  an  "electrical  index"  associated  with  the  lightning 
discharge  if,  in  fact,  they  do  not  represent  the  distribution  of  lightning  discharges 
over  the  Eastern  Hemisphere  for  1972. 

Heretofore,  the  lightning  discharge  distribution  has  been  assumed  to  follow 

3 

that  of  thunderstorm  days.  A comparison  of  the  estimates  of  these  distributions 
. produced  herein  for  January,  April,  July,  and  October  is  as  follows: 


I 

I 

i 

1 
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(1)  January,  The  correspondence  is  close  between  the  isoline  pattern  for 
January  of  thunderstorm  days  shown  In  Figure  4 and  that  of  lightning  discharges 
shown  in  Figure  13.  A broad  pattern,  centered  over  Madagascar  in  Figure  4 and 
South  Africa  in  P'igure  13,  extends  into  the  Indian  Ocean  in  both  charts.  Addition- 
ally, the  "low  center"  at  15N,  ISEin  Figure  4 is  present  in  Figure  13.  Other 
reasonably  close  correspondences  are  the  centers  in  both  charts  in  the  Mediter- 
ranean Sea,  Arabia,  India,  Japan,  and  Southeast  Asia.  The  patterns  in  Australia 
and  the  surrounding  oceanic  areas  differ.  (Again,  one  should  view  with  caution  the 
values  at  and  south  of  45S  in  all  charts  of  the  lightning  discharge  distribution.  ) 

(2)  April,  The  isoline  pattern  for  April  of  thunderstorm  days  shown  in 

Figure  6 and  that  for  the  number  of  lightning  discharges  shown  in  Figure  16  are 

dissimilar.  The  estimate  of  ten  thunderstorm  days  centered  at  35S,  45E  in 

Figure  6 is  associated  with  an  estimated  185,  000  lightning  discharges  in  F’igure  16 

2 

in  the  same  block  (1,  Oil,  540  km  ),  whereas  the  estimate  of  29  thunderstorm  days 
at  25N,  95E  in  Figure  6 is  associated  with  an  estimated  124,  000  lightning  discharges 
(1,  119,  165  km'^l  in  Figure  16.  As  explained  earlier,  this  type  of  reverse  corre- 
spondence caused  the  isoline  patterns  to  differ.  The  maximum  centers  of  30 
thunderstorm  days  in  Figure  6 at  15N,  45E  and  25N,  105E  have  somewhat  associ- 

4 

ated,  closed,  isolines  of  lightning  discharge  in  Figure  16  of  the  order  of  10  , but 
the  maximum  center  of  thunderstorm  days  at  5N,  85E  in  Figure  6 does  not  corre- 
spond to  a closed  pattern  of  lightning  discharge  in  Figure  16. 

(3)  July.  The  comparison  for  July  produced  mixed  results.  Although  both 
the  thunderstorm -day  results  shown  in  F'igure  8 and  the  lightning  discharge  results 
shown  in  Figure  19  have  centers  of  activity  over  Eurasia,  the  correspondence  of 
specific  centers  is  not  very  close.  The  thunderstorm -day  center  of  15  in  Figure  8 
at  15S,  115E  is  reflected  as  a lightning  discharge  center  in  Figure  19  of  the  order 

3 

of  10  to  the  northeast  of  Australia.  There  is  a very  close  correspondence  in  the 

thunderstorm  day  isoline  of  one  in  Figure  8 and  the  respective  lightning  discharge 
2 

isoline  of  10  shown  in  Figure  19,  both  of  which  span  the  northern  extremity  of 
the  Eurasian  continent. 

(4)  October.  There  is  a reasonably  close  correspondence  in  the  isoline  pat- 
terns of  estimated  thunderstorm  days  shown  in  Figure  10  and  that  of  lightning 
discharges  shown  in  Figure  22.  Both  patterns  have  high  centers  of  activity  over 
much  of  Africa,  the  Mediterranean  Sea,  and  India.  If  one  tracks  the  extensive 

Isoline  of  one  in  the  thunderstorm-day  chart  shown  in  Figure  10  and  compares  it 
2 

to  the  isoline  of  10  in  the  lightning  discharge  chart  shown  in  Figure  22,  it  becomes 
clear  that  the  respective  isoline  patterns  are  similar. 

Other  observations  concerning  the  estimated,  monthly  distributions  of  light- 
ning discharge  are; 
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(1)  The  center  of  tlie  estimated  number  of  lightning  discharges  in  South  and 
Central  Africa  is  either  the  dominant  center  (for  example,  in  January,  February, 
September,  October,  November,  and  December  shown,  respectively,  in  Figures  13, 
14,  and  21  through  24)  or  is  at  least  as  large  as  that  elsewhe-e  for  each  rntmth. 

(2)  In  all  the  months  that  the  combination  of  the  minimum  limits  of  threshold 
intensity  and  azimuthal  sectoring  permits  a judgment  (January  through  July  in 
Figures  13  through  10),  there  always  was  present  a center  of  lightning  discharge 
to  the  east  of  Australia  in  the  general  vicinity  of  the  Fiji  Islands.  There  also  is 
present  always  cither  a center  or  relatively  high  values  in  extreme  western 
Australia  and  the  coastal  waters  to  the  west. 

(3)  The  distribution  of  lightning  discharge  in  Southeast  Asia  is  not  as  sharply 
defined  in  terms  of  closed  isolines,  nor  is  it  as  intense  in  order  of  magnitude  as 
that  in  South  and  Central  Africa.  However,  there  always  is  present  a value  of  at 
least  10'’  lightning  discharges  in  Southeast  Asia. 

(4)  There  is  either  a center  of  lightning  discharge  or  an  isoline  of  at  least 

3 

10'  discharges  in  the  Mediterranean  Sea  for  each  month. 

(5)  Surprisingly,  there  is  a center  o',  lightning  discharge  or  an  Isoline  of  at 

3 

least  10  discharges  in  Arabia  for  each  month  except  July.  The  center  in  July, 

■> 

shown  in  Figure  19,  is  only  of  the  order  of  10'  discharges. 

(f!)  One  may  follow  the  effect  of  the  march  of  the  seasons  in  the  Northern 

Hemisphere  (0-179E)  by  looking  at  the  general  latitudinal  position  of  the  isoline  of 
2 

10  lightning  discharges  that  spans  Furasia  in  each  month,  as  shown  in  Figures  It 
through  24.  In  general  terms,  there  is  a march  poleward  of  the  10^  isoline  from 
December  to  July.  In  July,  as  shown  in  Figure  lit,  there  are  two  values  in  the 
Arctic  at  85N  latitude.  The  isoline  of  10‘  discharges  gradually  recedes  equator- 
ward  from  July  to  the  end  of  fall  in  December. 

The  estimated  number  of  lightning  discharges  for  1972  is  shown  in  Figure  25. 
The  broken  line  is  the  threshold-sector  line  for  the  period  January  to  July  and  the 
dash-dot  line  is  the  threshold-sector  line  for  the  period  August  to  December. 
Therefore,  one  should  realize  that  only  the  values  north  of  the  dash-dot  line  were 
based  on  data  for  all  of  1972.  It  is  to  be  recalled  that  values  south  of  4r).S  were 
questionable  since  they  were  based  on  only  a relatively  few,  intense,  recorded 
signals  and  the  assumption  that  the  log-normal  probability  distribution  was  valid 
throughout  the  range  of  possible  field  strengths,  including  the  extreme  tails  of  the 

4 

distribution.  It  seems  unlikely  that  the  10  isoline  of  lightning  discharge  in  the 
•Southern  Ocean  and  over  Antarctica  is  valid. 

The  focus  of  thunderstorm  days  in  Africa  also  is  reflected  in  the  distribution 
of  lightniiig  discharges  shown  in  Figure  25.  The  highest  estimate  of  1,  2H3,  300 
discharges  over  some  1,  230,  1(>3  km^  for  all  of  1972  is  located  at  5S,  15K  in  the 

3 

focus  of  South  Africa.  There  is  a center  of  lightning  discharge  of  the  order  of  10' 
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(lischarge.s  for  the  year  over  the  Me(Jiterranean  Sea.  The  foeu-s  of  thunderstorm 

days  in  Southeast  Asia  was  weaker  than  that  of  Africa  both  in  this  study  and  in  the 

\VM()  compilation.  The  estimated  number  of  lightning  dischai  ges  in  Southeast 

Asia  likewise  was  weaker  than  the  African  estimate  by  an  order  of  magnitude  in 
4 5 

general  (that  is,  10  vs  10  ) and  by  two  orders  of  magnitude  for  the  respective 

central  values  (1,  283,  300  at  fiS,  15K  vs  42,  400  at  15N,  OfiK).  There  were  relatively 

4 

large  estimates  for  the  number  of  lightning  discharges  (of  the  order  of  10  ) in 
Western  Australia  and  over  the  islands  in  the  South  Pacific  to  the  northeast  of 
Australia. 

If  one  accepts  only  values  north  of  43S  latitude,  as  suggested  earlier,  it  is 

4 

interesting  to  compare  the  areal  extent  mapped  out  by  the  isoline  of  10  discharges 
in  I-'igure  25  with  the  region  in  the  Kastern  Hemisphere  which  is  influenced  by  the 

monsoons.  As  mentioned  earlier,  Hamage'  *’  delineates  the  monsoon  region  as  that 

4 

enclosed  between  3.5.\  and  25S  and  between  30\V  and  170K.  The  isoline  of  10  dis- 
charges north  of  45S  in  Figure  25  and  the  monsoon  region  after  Hamage  for  the 

Kastern  Hemisphere  are  reproduced  in  Figure  20.  As  we  can  see  from  Figures  25 

4 3 f) 

and  2fi,  the  isolines  of  the  highest  three  orders  of  magnitude  (10  , 10  , and  10  ) 
of  the  estimated  number  of  lightning  discharges  were  enclosed  within  the  monsoon 
region  to  a significant  degree.  This  was  not  surprising  since  the  tropics  are  con- 

3 

sidered  the  dominant  source  of  thunderstorms  on  Earth.' 

!1.  1 111*'  Int  nf  l.i^hlnint:  lli-t-hurgr 

A few  investigators  have  given  point  or  areal  estimates  of  the  incidence  of 
lightning  discharge.  These  estimates  are  presented  in  Table  9.  As  used  here, 
the  incidence  of  lightning  discharge  is  defined  as  the  occurrence  of  lightning  dis- 
charges during  a stormy  period  per  unit  area  per  unit  time.  The  data  of  this  study 
were  amenable  to  an  analysis  of  the  distribution  of  the  incidence  of  discharge  for 
much  of  the  Eastern  Hemisphere. 

As  in  most  studies  in  meteorology,  certain  assumptions  must  be  made.  Fol- 
lowing Horner,  **  the  average  extent  of  a "thunderstorm  area"  was  assumed  to  be 
1000  km^;  that  is,  during  periods  when  thunderstorms  were  active  in  a data  block, 
the  average  areal  extent  of  the  storm  (or  storms)  was  assumed  to  be  1000  km‘. 

The  duration  of  thunderstorms  on  the  average  was  assumed  to  be  3 hr  in  the 
tropics**^  and  I hr  in  temperate  and  polar  latitudes.^  The  "tropics"  was  somewhat 
arbitrarily  defined  as  encompassing  the  area  between  20N  and  20.S. 

The  incidence  of  lightning  discharge  at  each  grid  point  was  estimated  for  the 
classic,  climatologic  months  of  January  and  April.  First,  an  estimate  of  the 
number  of  discharges  per  thunderstorm  day  for  the  month  was  obtained.  Each 
grid  point  in  the  distribution  of  the  estimated  number  of  lightning  discharges  for 
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Figure  2(i.  A Comparison  of  the  Area  Knclosed  within  the  Isolines  of  10 
Lightning  Discharges  in  the  Kastern  Hemisphere  for  1972  North  of  45.S 
I.atitude  (heavy  solid  line)  and  the  Hegion  that  is  I’sually  under  the  Influ- 
ence of  the  Monsoons  (shad<-d  box) 


Table  9.  Kstmiates  of  Hie  Incidence  of  Lightning  Discharge  from  the  Literature 
(Lnits  of  km"^  sec"*) 


Incidence  of 
Lightning  Discharge 
(km"*^  sec*’) 


Source 


Alva  (1968) 


Main  thunderstorm  areas 


Horner  (1965) 


Hrooks  (1925) 


The  main  thunderstorm  areas  are  the  Americas,  Africa,  and  Southeast  Asia 
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January  and  April,  shown  in  Figures  13  and  16,  respectively,  was  divided  by  the 

corresponding  estimate  of  the  number  of  thunderstorm  days  shown  in  Figures  4 

and  6.  This  "resultant  matrix"  was  divided  by  either  10,  800,  000  for  the  tropics 

(that  is,  1000  km  X 10,  800  sec)  or  3,  600,  000  elsewhere  (that  is,  1000  km  X 

3,  600  sec).  The  final  result  at  each  grid  point,  shown  in  Figure  27  for  January 

and  Figure  28  for  April,  was  an  estimate  of  the  incidence  of  lightning  discharge  in 
-2  -1 

units  of  km  sec 

Again,  values  at  and  south  of  45S  were  questionable.  The  values  in  Figures  27 
and  28  are  negative,  ordinary  logarithms  of  the  incidence  of  lightning  discharge. 

Of  Course,  the  mantissa  for  each  value  is  one  less  than  the  appropriate  negative 
power  of  10.  For  example,  in  Figure  27  at  15N,  125E,  one  finds  the  value  4.81 

and  may  immediately  surmise  that  the  estimate  of  incidence  for  this  data  block  is 
-5  -2  -1 

of  the  order  of  10  lightning  discharges  km  sec  during  a stormy  period. 
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Figure  27.  Distribution  of  the  Incidence  of  Lightning  Discharge  (km'^  sec'S  for 
January  1972.  The  dash-dot  line  is  the  threshold-sector  line  for  the  period  Janu- 
ary to  July  1972.  Values  are  negative,  ordinary  logarithms  for  10-deg  data  blocks 
centered  at  the^  points  shown. 
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Figure  28.  Distribution  of  the  Incidence  of  Lightning  Discharge  (km  sec  ) for 
April  1972,  The  dash-dot  line  is  the  threshold-sector  line  for  the  period  January 
to  July  1972.  Values  are  negative,  ordinary  logarithms  for  10-deg  data  blocks 
centered  at  the  points  shown 

Estimates  of  the  incidence  of  lightning  discharge  by  latitude  zones  are  given  in 
Table  10  for  January  1972  and  Table  11  for  April  1972.  The  respective  estimates 

for  the  two  months  agree  consistently  in  order  of  magnitude,  although  they  vary 

-5  -2  -1 

from  a close  agreement  of  96  percent  (that  is,  2,3  X 10  km  sec  compared  to 
-5  -2  -1 

2.  4 X 10  km  sec  ) for  the  estimates  in  the  zone  of  high  north  latitudes  to  a 

disagreement  by  a factor  of  3 for  the  estimate  of  all  available  data  for  each  month 

(that  is,  3.  6 X 10  ^ compared  to  1.  1 X 10'^  km”^  sec*^). 

Aiya'^  estimated  that,  in  general,  the  incidence  of  lightning  discharge  in  India 
-5  -2  -1 

is  9 X 10  km  sec  as  indicated  in  Table  9.  The  estimates  on  the  average  of 

-5 

this  study  for  those  data  blocks  that  include  some  portion  of  India  were  3.5X10 
-2  -1 

km  sec  for  both  January  and  April  of  1972,  The  estimates  of  this  study  for 
that  portion  of  the  tropics  in  the  Northern  and  Eastern  Hemisphere,  (that  is,  be- 
tween the  equator  and  20N  and  between  the  Greenwich  Meridian  and  179E)  were 
1.6X10  ^ km'^  sec  ^ for  January  1972  and  3.  4 X lO"^  km'^  sec*^  for  April  1972, 
as  shown  in  Tables  10  and  11,  respectively. 
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Table  10.  Estimate  of  the  Incidence  of  Lightning  Discharge  by  latitude  Zone  for 
January  1972  (Units  of  km"2  sec'^) 


Latitude  Zone 

Number  of 
Values 

Sum  of  Values 
(km'2  sec"*) 

Average  Value 
(10"5  km"2  sec")) 

High  north  latitudes 
(60N-90N) 

3 

6.  9 X 10"^ 

2.3 

Middle  latitudes 
(21N-59N) 

48 

3.  0 X 10'^ 

6.3 

N.  H.  tropics 
(0-20N) 

34 

5.  5 X 10“* 

1.6 

S.  H.  tropics 
(0-203) 

34 

2.  8 X 10"^ 

8.3 

(Mher 

(part  of  21S-9LS, 
see  Figure  27) 

72 

2. 8 X 10"^ 

38 

N.  Hemisphere 
portion  of 

E.  Hemisphere 

85 

3.4  X 10"^ 

4.  2 

Eastern  Hemisphere 
tropics  (20N-20S) 

68 

3. 4 X 10'^ 

5.0 

All  grid  point 
estimates* 

191 

3. 4 X 10'^ 

36 

* 


Includes  values  along  and  south  of  45S  latitude, 


which  may  be  questionable. 


Horner  gave  an  order  of  magnitude  estimate  of  the  incidence  of  discharge  of 
-5  -2  -1 

10  km  sec  for  the  "main  thunderstorm  areas  of  the  world'  as  shown  in 
Table  9.  A roughly  comparable  estimate  is  that  of  this  study  for  the  incidence  of 
discharge  in  the  tropics  of  the  eastern  Hemisphere  (that  is,  20N  to  20S)  as  given 
in  Table  10  for  January  1972  and  Table  11  for  April  1972.  These  estimates  were 
5.0  X 10  ^ km  ^ sec  ' for  January  and  2.  (i  X 10  ^ km  ^ sec  * for  April. 

Horner  also  gave  a range  of  estimates,  based  on  various  techniques  men- 
tioned in  Section  1,  for  the  global  incidence  of  lightning  discharge.  These  estimates 
•6  *2  *1  "6  *2  *1 

ranged  from  6X10  km  sec  ‘tol6X10  ” km  ^ sec  * as  shown  in  Table  9. 

One  may  infer  from  Brook's  work,  as  mentioned  in  Section  1 and  presented  in 

-4 

Table  9,  that  the  global  incidence  of  lighining  discharge  is  of  the  order  of  10 
-2  -1 

km  sec  . The  most  reliable  planetary-scale  estimates  of  this  study  of  the 

incidence  of  discharge  were  for  the  Northern  Hemisphere  (0-179E).  The  estimate 

-5  -2  -1 

of  the  incidence  for  January  1972  for  this  area  was  4.  2 X 10  km  sec  as  given 
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TabU’  11.  Kstiniato  (if  the  Incuience  of  l.ightning  Discharge  by  latitude  '/.onp  for 
April  1P72  (I  nits  of  km'-  sec'*) 


l .atituile  Zone 

Number  of 
Values 

■Sum  of  Values 
(km'^  sec'l) 

Average  Value 
( lO"-’  km'2  sec'  ^ ) 

High  north  latitudes 
(r>0N-90N) 

4 

9.  (!  X lO'''’ 

2.  4 

Middle  latitudes 
(21N-.i9N) 

.39 

1.  7 X 10*'^ 

2.  8 

H.  tropics 
(0-20N) 

3 5 

1. 2 X lO"'^ 

3.4 

■S.  H.  tropics 
(0-20.S) 

3(1 

1.  5 X 10"^ 

4.  2 

( ither 

(part  of  21.S-91.S, 
sec  Figure  2H) 

70 

1. 8 X 10"^ 

2.  5 

N.  Hemisphere 
portion  of 

K.  Hemisphere 

OH 

2.  9 X lO"'^ 

3.  0 

Kastern  llemispliere 
tropics  (20N-20.S) 

10(1 

2.  7 X lO'* 

2.  r, 

All  grid  jx>int 
estimates* 

204 

2. 2 X 10‘^ 

1. 1 

* 


Includes  values  along  and  south  of  4 5S  latitude. 


which  may  be  questionable. 


in  Table  10  and  .3.0  X 10  km  ^ sec  ^ for  April  1972  as  given  in  Table  11.  The 

estimate  of  the  average  incidence  of  discharge  for  all  grid  point.s  that  were  avail- 

-4  -2  - 1 

able  above  the  threshold-sector  line  in  January  1972  was  3.  fi  X 10  km  sec  and 

-4  -2  -1 

the  commensurate  estimate  for  April  1972  was  1.  1 X 10  km  sec  as  indicated 
in  Tables  10  and  11,  re.spectively.  However,  values  in  high  latitudes  of  the 
■Southern  Hemisphere  that  were  used  in  these  estimates  were  questionable. 

3. .7  Nri'iil  roni'cntriilion  (»f  l.lahlning  llisrhargi* 

There  seems  to  be  no  e.stlmate  available  in  the  literature  for  the  areal  concen- 
tration of  lightning  discharge.  This  parameter  of  thunderstorm  activity,  as  used 
here,  is  the  number  of  lightning  discharges  per  unit  area  over  some  time  period 
irrespective  of  the  specific  areal  extent  of  thunderstorm  activity  which  caused  the 
discharges.  In  other  words,  it  is  simply  the  number  of  discharges  per  gross  area. 


'Ihis  paraniotcr  is  no^  (hr  .sanir,  thrrrforr,  as  (hr  "lightning  discharge  density," 
which  is  always  related  to  the  extent  of  area  iinrier  the  influence  of  thunderstorms. 

An  estimate  of  the  distribution  of  the  areal  concentration  of  lightning  discharge 
for  1972  was  obtained  by  divdding  the  value  at  each  grid  point  in  the  yearly  estimate 
of  the  total  number  of  lightning  discharges  (f  igure  25)  by  the  area  of  the  10-deg, 
latitude-longitude  block  that  encloses  the  grid  [Joint.  The  areas  for  the  nine  blocks 
from  pole  to  eqtiator  in  a hemisphere  arc  listed  in  Ttible  12.  These  values  were 
calculated  by  a straightforward  application  of  spherical  trigonometry  on  the  basis 
of  a spherical  Karth  (r  (>371  km). 


Taljle  12.  Areas  on  Karth  Hounded  by  Any  Two  Meridians  10-deg  Apart  and  the 
Two,  Indicated  Parallels;  the  Average  Areal  Concentration  of  Lightning  Discharge 
in  the  Northern  Hemisphere  Within  the  Indicated  Parallel  Hounds 


A rea 

Num  ber 

, 

Parallel  Hounds 
(degrees) 

Area 

(km^) 

Average  Areal 
Concentration 
of  Discharge  (km'^) 

1 

0-10 

1,  230,  Ui3 

1.31  X 10'* 

2 

10-20 

1,  192,  786 

3.  94  X 10'^ 

3 

20-30 

1,  1 19.  165 

3.48  X 10'“ 

4 

30-40 

1 . 001.  540 

4.09  X lO'^ 

5 

0 

1 

6 

873,  180 

1. 27  X 10'‘* 

6 

50-60 

708,  288 

6.  73  X lO*'^ 

7 

0 

1 

O 

521,  875 

2.  68  X 10"'^ 

a 

70-80 

310,  606 

8.  58  X 10“* 

-4 

9 

80-90 

107,  625 

9. 33  X 10 

The  yearly  estimate  of  the  distribution  of  the  areal  concentration  of  lightning 
discharge  is  given  in  figure  29.  All  values  are  negative,  ordinary  logarithms  of 
the  areal  concentration  of  discharge  except  the  value  rjf  0.02  at  5S,  15K  which  is 
positive  (that  is,  the  logarithm  of  1,  2K3,  300  discharges  at  5S,  IfjK  from  figure  25 
divided  by  1.  230,  163  km^  from  Table  12.)  T'he  threshold-sector  line  and  the  con- 
comitant explanation  are  the  same  for  the  areal  concentration  of  discharge  as  that 
for  the  yearly  estimate  for  1972  of  the  number  of  lightning  discharges  given  in 
figure  25  and  on  page  59.  Kecall  (hat  only  the  values  in  figure  25  north  of  the 
dash-dot  line  are  based  on  data  for  all  of  1972. 
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rigure  2!>.  Distribution  of  tiic  Areal  Concentration  of  Lightning  Discharge  (km  ) 
for  1972.  The  broken  line  is  the  average  threshold-sector  line  for  the  period 
January  to  July  1972  and  the  dash-dot  line  that  for  the  period  August  to  December 
1972.  Values  are  negative,  ordinary  logarithms  (except  for  that  at  .'iS,  1 r>K  which 
U positive)  for  lO-deg  data  blocks  centered  at  the  points  shown.  Solid  box  is  the 
area  of  maximum  data  coverage  in  the  Northern  llemis[)here 

One  may  approximate  the  latitudinal  change  of  the  ai-eal  concentration  of 
lightning  discharge  from  a study  of  the  area  of  maximum  data  coverage  in  the 
Northern  Ilemiisphere  (I'igure  2!t).  as  outlined  by  the  solid  box  (that  is,  by  the 
grid  points  of  5N  to  !(r)N  latitude  and  .‘jL  to  lOSlC  longitude).  The  average  areal 
concentration  of  discharge  for  the  10-deg  zones  of  latitude  with  central  points  as 
shown  in  I'igure  2!'  (that  is,  fr.N.  Ki\,  2r)N,  and  so  on)  is  presented  in  Table  12. 
The  general  trend  r-vident  from  Table  12  and  I'igure  29  was  that  the  average,  areal 
conrarnt  rat  ion  of  lightning  dischargr;  'U'creased  with  geographic  latitude  in  1972  in 
the  .Northern  Hemiipherc  quadrant  of  the  Lastern  Hemisphere.  The  areal  concen- 
tration of  discharge  wa-.  three  orders  of  magnitude  higher  on  the  avio-age  in  the 
equatorial  zone  (O-ION)  than  in  the  |K)lar  zone  (H0\-90.N)  for  ltt72. 


S.t>  Niun  r-  »f  ^.r^llr 


There  were  various  possible  sources  of  error  in  the  initial  c(jllection  of  data 
and  subsequent  analyses,  some  of  which  have  been  mentioned  previously.  The 
known,  possible  sources  of  error  were: 

(1)  The  instrument  error  is  not  known  in  quantitative  terms.  However,  errors 
in  calibration,  the  setting  of  the  gating  circuits,  and  so  on,  are  possible.  More- 
over, two  signals  that  are  incident  at  the  antenna  at  precisely  the  same  instant  in 
time  would  cause  an  erroneous,  distorted  waveform  to  be  recorded.  The  time 

-5 

resolution  of  the  signal  monitor  for  individual,  received  sferics  is  10  sec. 

(2)  The  threshold  arrays  (lower  limits  on  the  probability  integrals)  rre  only 

.32 

mean  estimates,  dependent  to  a large  extent  on  propagation  laws  at  VI. F.  Fierce 
stated  that  propagation  at  VI. F is  subject  to  geographic,  geomagnetic,  and  tem- 
poral variability,  none  of  which  is  known  with  certainty. 

(3)  The  important  statistics  listed  in  Table  2,  which  generally  are  point  esti- 
mates in  the  mean,  may  not  be  suitable  for  this  work.  Since  there  is  wide  agree- 
ment that  the  standard  deviation  of  received  sferics  ranges  from  (i  to  8 dB,  we  may 
appreciate  the  effect  on  the  analyses  of  varying  this  sample  statistic,  other  things 
being  equal. 

An  estimate  of  the  distribution  of  the  number  of  lightning  discharges  for 
August  1972  with  standard  deviations  of  (1,  7,  and  8 dH  are  shown  in  h’igure  30. 

The  middle  distribution  in  Figure  30,  presented  in  hundreds  of  discharges  with  a 
standard  deviation  of  7 dll,  is  the  same  as  that  given  previously  in  Figure  20  for 
August  1972.  For  a given  threshold  value  at  a grid  poin^  (that  is,  a given  lower 
limit  on  the  probability  integrals  as  shown  in  Kq.  (13)1,  the  estimate  of  the  number 
of  lightning  discharges  increases  with  increasing  standard  deviation  as  shown  in 
Figure  30  (or  remains  the  same  duo  to  round-off  to  the  nearest  100  discharges). 
Since  standard  deviation  is  a measure  of  dispersion,  the  probability  of  occurrence 
above  a given  threshold  value  (namely,  ftercentile)  increases  with  increasing 
standard  deviation.  Hut  the  probability  integrals  that  are  thus  effected  enter  into 
Fq.  (13)  in  the  denominator  with  the  result  that  the  estimate  of  the  number  of 
lightning  discharges  decreases  at  a grid  point  with  an  increase  in  standard 
deviation. 

A detailed  analysis  of  error  of  the  values  in  Figure  30  is  not  warranted  be- 
cause of  the  uncertainties  associated  with  the  other  sample  statistics  given  in 
Table  2 that  affect  the  distribution  of  sferics.  In  general,  there  is  little  or  no 
change  in  the  estimates  of  the  number  of  lightning  dischai'gcs  (to  the  nearest  100 
discharges)  in  high  north  and  south  latitudes  wliere  the  threshold  array  places  the 
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Fiifure  30.  Comparison  of  the  Distributions  of  l.i^htnitiR  Disciiargo  for  August  1972 
for  Standard  Deviations  of  Ucceived  Sferics  of  (a)  (i  dK,  (b)  7 dll,  and  (c)  8 dI3. 
Values  are  in  hundreds  of  discharges  for  lO-deg  data  blocks  centereil  at  the  points 
shtjwn.  The  broken  line  is  the  average  threshold-sector  line  for  the  period  August 
to  December  1972 
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appropriate  t value*  in  the  extreme  "left"  and  "right"  tails  of  the  respective  dis- 
tributions. The  greatest  change  in  results  caused  by  varying  the  standard  devia- 
tion alone  is  roughly  50  percent.  For  example,  at  5N,  5E  the  estimates  of  32,  400 
discharges  for  a standard  deviation  (SD)  of  6 dB  is  47  percent  greater  than  the 
estimate  of  15,  100  discharges  for  a SD  of  7 dB;  and  the  estimate  of  8600  discharges 
for  a SD  of  8 dB  is  57  percent  less  than  that  for  a SD  of  7 dB. 

(4)  The  distributions  of  the  various  types  of  sferics  studied  may  not  be  log- 
normal over  certain  ranges  of  values  of  the  respective  random  variables.  Indeed, 
these  sferics  may  follow  some  other  distribution  for  extensive  periods  of  time. 

(5)  The  analyses  of  the  number  of  discharges,  the  incidence  of  discharge,  and 
the  areal  concentration  of  discharge  depend  for  efficacy  on  the  rough,  empirical 
equations  quoted  or  developed  in  this  section.  Moreover,  these  analyses  depend 
on  assumptions  that  may  not  be  valid. 

The  various  possible  sources  of  error  mentioned  above  were  very  likely  of  a 
random  nature  and  may  have  compensated  for  each  other  to  some  degree. 

».  CONCl.l  SHINS  \M)  HK«  tlMMKNDXTIONS 

l.l  l!onc'lusion> 

The  following  distributions  of  parameters  were  developed  in  this  study: 

(1)  thunderstorm  days  for  January,  April,  July,  and  October  of  1972  and  for  the 
entire  year;  (2)  the  number  of  lightning  discharges  for  each  month  of  1972  and  for 
the  entire  year;  (3)  the  incidence  of  lightning  discharge  for  January  and  April  of 
1972;  and  (4)  the  areal  concentration  of  lightning  discharge  for  1972,  The  study 
was  restricted  to  the  Eastern  Hemisphere  and  the  data  were  unique  in  that  no 
sferics  data  on  a planetary-scale  basis  have  been  heretofore  available.  Therefore, 
the  analyses  based  on  these  data  resulted  in  first  estimates  of  the  distribution  of 
the  above  parameters,  except  for  that  of  thunderstorm  days.  There  were  few 
uncertainties  associated  with  the  estimates  of  thunderstorm  days,  but  there  were 
many  uncertainties  in  the  other  analyses  and,  therefore,  the  results  of  these  anal- 
yses must  be  considered  tentative.  Nevertheless,  the  following  conclusions  for 
the  Eastern  Hemisphere  in  1972  may  be  stated; 

(1)  The  thunderstorm-day  charts  of  this  study  differ  considerably  from  the 
mean,  monthly  charts  of  the  WMO  that  were  based  on  data  compiled  in  1953.  One 
possible  cause  of  the  differences  was  that  the  monthly  charts  of  this  study,  that 


Where  t • (Tj,  - T)/a;  Tq  and  T are  the  threahold  intensity  in  field  strength  at  a 
grid  point  and  the  mean  field  strength  of  the  distribution,  respectively;  o is 
standard  deviation. 
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wero  based  on  one  year  of  data,  could  represent  an  anfinialy  from  the  mean  pattern 
of  the  WMO  charts.  Another  cause  of  difference  that  seems  incontrovertible  was 
that  in  the  WMO  compilation  of  the  data  there  were  extensive  areas  on  Karth  with 
a void  of  data  coveraf^e.  T hese  areas  include  the  oceans.  North  Africa,  Arabia, 
the  Mediterranean  Sea,  and  Southeast  China,  ft  has  been  show'n  in  this  study  that 
there  were  many  more  thunderstorm  days  in  these  areas  in  1972  than  reported  by 
the  WMO  in  the  mean. 

(2)  A first  estimate  of  the  distribution  of  the  number  of  lightning  discharges 
over  a large  area  of  Karth  was  presented  in  Figures  13  through  25.  The  validity 

of  the  magnitude  of  any  one  estimate  at  a grid  point  or,  indeed,  of  all  the  estimates, 
is  not  known.  However,  the  magnitude  of  the  values  relative  to  each  other  and  the 
isoline  patterns  probably  represent  a useful  index  of  "electrical  activity  over  a 
broad  area  of  Karth. 

(3)  There  was  a close  correspondence  between  the  isoline  patterns  of  thunder- 
storm days  and  lightning  discharges  for  the  months  of  January  and  October.  There 
was  little  correspondence,  in  general,  between  these  isoline  patterns  for  April  and 
July.  (Heretofore,  the  assumption  has  been  made  frequently  that  the  distributions 
w ere  the  same. ) 

(4)  Central  and  South  Africa  were  the  dominant  foci  of  lightning  occurrence  on 
a yearly  basis  and  either  a relatively  high  center,  or  the  highest  center  of  activity 
on  a monthly  basis.  The  occurr  ence  of  lightning  was  relatively  high  in  the 
Mediterranean  Sea  and  Arabia  for  most  months.  The  occurrence  of  lightning  in 
Southeast  Asia  was  consistently  smaller  in  magnitude  than  in  Africa.  There  was 
always  present  a relatively  high  occurrence  of  lightning  in  Kastern  and  Uestern 
Australia  and  in  the  adjacent  oceanic  areas. 

(5)  The  latitudinal  occurrence  of  lightning  discharges  was  shown  to  follow  Sun 
northward  over  Kurasia  from  January  to  a maximum  poleward  thrust  in  July  with 
occurrences  of  discharge  over  the  Arctic.  The  occurrence  of  discharge  then 
receded  equatorward  from  the  Arctic  from  July  to  the  end  of  fall  in  December. 

(ti)  An  estimate  of  the  distribution  of  the  incidence  of  lightning  discharge  for 
January  and  April  was  presented  in  bigures  27  and  28,  respectively.  A breakdown 
of  the  average  incidence  of  discharge  by  zones  of  latitude  was  presented  in 
'l  ables  10  and  11.  1 he  most  reliable  planetary-scale  estimate  of  the  average 

incidence  of  discharge  was  4.2  X 10  ^ km  ^ sec  ^ for  January  and  3.0  X 10  ^ km 
sec  ^ for  April  over  the  Northern  Hemisphere  (0-179K).  An  average  incidence  of 
discharge  of  3.  (i  X lo“*  km'^  sec*^  was  obtained  for  January  and  1.1X10  ^ sec  * 
for  April  over  the  Kastern  Hemisphere  as  a whole,  but  these  results  were  con- 
sidered less  reliable  due  to  uncertainties  in  the  analysis. 

(7)  An  estimate  of  the  distribution  of  the  areal  concentration  of  lightning  dis- 
charge was  presented  in  Figure  29.  There  was  evidence,  based  on  an  analysis  in 
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the  Northern  Hemisphere  (0-105E),  that,  in  general,  the  areal  concentration  of 
discharge  decreased  with  geographic  latitude. 

I.s  Ki-i'ommi-ndiilion- 

(1)  There  is  much  current  interest  in  the  relationship  between  thunderstorm 

32 

days  and  the  incidence  of  lightning  discharge.'  A direct  extension  of  this  work 
would  be  the  development  of  laws  that  relate  these  parameters  on  a monthly  and 
yearly  basis. 

(2)  Very  useful  information  on  the  diurnal  variation  of  thunderstorm  activity 
could  be  developed  from  these  data.  This  is  especially  true  for  oceanic  areas 
where  frequently  the  diurnal  variation  is  assumed  to  be  constant.  Kstimates  of 
the  diurnal  variation  of  thunderstorm  activity  for  Africa,  Kurasia,  Southeast  Asia, 
Australia,  and  the  Kastern  Hemisphere  as  a whole  could  be  developed. 

(3)  One  could  analyze  the  extensive  data  in  (>-hr  increments  to  obtain  estimates 
of  the  number  of  active  thunderstorms  that  occur  on  Karth  at  any  one  time.  This 
information  may  be  useful  to  those  who  study  the  general  circulation  of  the 
atmosphere. 

(4)  An  analysis  of  the  data  could  be  performed  in  (i-hr  increments  and  esti- 
mates made  of  the  disturbed-weather  electrical  field  over  a large  area  of  Karth 
that  is  generated  by  thunderstorms. 

(5)  One  could  use  both  visible  and  infrared  satellite  photographs  for  1972  that 
show'  discernible,  large,  isolated,  cumulonimbus  clouds  and  a fi-hr  "electrical 
index"  developed  from  these  data  to  attempt  to  relate  the  sferics  count  to  the  ver- 
tical and  areal  extent  of  the  cloud  shield.  A strong  correlation  could  lead  to  a 
useful  analysis  of  the  electrical  intensity  of  thunderstorms  as  it  relates  to  the 
diameter-to-depth  ratio  of  clouds. 

(fi)  A challenging  investigation  would  be  to  attempt  to  fit  an  empirical,  statis- 
tical distribution  to  the  occurrence  of  lightning  discharge  on  a fi-hr  basis  over  a 
large  area  of  Karth.  These  data  are  amenable  to  such  investigation. 

(7)  (me  could  build  evidence  based  on  these  data  directed  at  showing  that  nature 
dominantly  selects  the  ice  process  in  cloud  electrification  over  that  of  liquid  water. 
The  steps  in  this  investigation  would  be  first  to  develop  a bench  mark  of  electrical 
activity,  K,  over  a grid.  This  could  be  done  by  normalizing  the  sferic  counts  by 
the  use  of  the  threshold  arrays.  Then  one  could  regress  E onto  latitude  and  per- 
haps obtain  a rough  fit  to  the  data.  Next,  one  could  find  a best -fit  approximation 
for  E in  terms  of  geographic  latitude  and  the  mean  terrain  height.  It  is  hoped  that, 
the  fit  would  be  better  than  previously,  finally,  one  could  find  the  best-fit  approx- 
imation of  E in  terms  of  latitude,  0,  mean  terrain  height,  H^,  and  the  freezing 
level,  FI,,  of  the  form 
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K r(0.  Hj.  FL)  . (ir,) 

If  one  could  aaserf,  with  a high  level  of  confidence,  that  this  best-fit  approximation 
was  a statistically  significant  improvement  over  the  previous  ones,  then  attention 
could  be  directed  to  developing  a physical  explanation  of  why  the  improvement 
came  about. 
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Remote  Sensing  Technique 


The  purpose  of  this  appendix  is  to  expand  the  information  given  in  Section  2 on 
the  sferics  network.  The  network  is  comprised  of  signal  monitors,  which  are 
electronic  systems  that  detect  and  measure  random  occurring  electromagnetic 
signals  (sferics)  produced  by  lightning  discharges.  Ten  signal  monitors  were  in 
operation  during  the  period  .lanuary  to  July  1972  and  six  during  the  remainder  of 
the  year.  Six  sensors  arc  sufficient  to  monitor  and  record  sferics  on  a global 
basis.  As  mentioned  previously,  this  study  was  confined  to  the  Kastern  Hemisphere 
The  signal  monitor  equipment  and  recording  methodology,  sferics -fix  techniques, 
and  the  system  threshold  arrays  will  be  explained.  Although  available,  details  of 
the  electronic  subsections  of  the  signal  monitors  are  considered  beyond  the  scope 
of  this  study  and  will  not  be  included. 


\l.  'I(,N  \l  MUM  KIH  l yi  l|•Ml:^i 

The  signal  monitor  ecpiipment  is  depicted  in  Figure  Al.  A vertical,  non- 
directional  whip  antenna  and  two  directional  loop  antennae  are  used  to  receive  the 
sferic  waveform.  Although  all  electromagnetic  signals  appear  on  the  antennae, 
the  signal  monitors  record  on  a selective  basis,  using  the  criteria  of  waveform 
half-cycle  characteristics,  the  relative  threshold  field  intensity  of  the  waveform, 
and  the  azimuth  of  arrival.  The  signal  monitor  also  has  the  capability  of  recording 
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Figure  Al.  Signal  Monitor  Operation 


liignal.-s  selectively  based  on  the  energies  of  the  distributed  frequencies  that  com- 
pose the  waveform.  However,  the  spectrum  analyzer  section  and  the  discriminant 
section  which  perform  this  function  were  not  in  operation  during  1972.  'I'he  individ- 
ual waveforms,  after  analysis  by  the  waveform  and  azimuth  sections,  were  recorded 
to  an  accuracy  of  10  ^ seconds. 

A typical  recorded  waveform  is  shown  in  Section  2 (Figure  1).  The  waveform 
has  the  appearance  of  a sine  wave  which  builds  to  a peak  and  then  decays  to  zero  in 
1000  unec  or  less.  The  half-cycle  of  maximum  amplitude  of  the  recorded  wave- 
form occurs  within  the  first  .'100  fjscc  of  the  waveform.  The  peak  energies  of  all 
recorded  sferics  occur  in  the  frequency  range  of  250  Hz  to  fiO  kHz. 

The  whip  antenna  has  a precise  electrical  length  (2  m)  for  establishing  the  rela- 
tive field  intensity  in  volts  per  meter  for  each  received  waveform,  irrespective  of 
direction  of  propagation.  Klectromagnetlc  signals  appearing  across  the  antenna 
produce  an  output  voltage  proportional  to  its  electrical  height  and  the  received  sig- 
nal intensity.  The  waveform  input  circuits  receive  the  output  of  the  vertical  whip 
antenna.  The  signal  monitor  input,  threshold  control,  peak  detector  and  gate 
selector,  and  normalizing  circuits  in  conjunction  with  the  loop  antennae  and  azi- 
muth sectirjn  circuits  perform  the  prerecording  criteria  determinations  previously 
described. 


HO 


Tho  loop  antennae  are  u.scd  to  resolve  the  received  si^^nal  into  two  waveforms 
having  an  angular  separation  of  90  degrees.  The  two  waveforms,  in  turn,  produce 
coordinate  voltages  that  are  used  to  determine  the  received  azimuth  of  the  signal 
in  degrees  relative  to -true  north.  The  azimuth  accuracy  approximates  one  degree 
with  a full  3(i0-degree  coverage. 


\1.  'I  I UK  I l\l\(,  li  t HMtJI  I 

The  source  location  on  earth  of  a sferic  has  been  determined  by  the  use  of 
recording  stations.  'I'hc  time  of  arrival  of  the  sferic  at  each  station  was  recorded 
to  an  accuracy  of  10  sec.  Them,  tho  difference  in  time  of  arrival  between  sta- 
tions was  Converted  to  distance  by  the  use  of  the  propagation  rate  at  VI, K.  Although 
empirical  tests  were  conducted  to  determine  an  eeptimum  propagation  rate,  it  was 
found  that  the  speed  of  light  was  acceptable.  The  speed  of  light  is  the  speed  of 
propagation  of  electromagnetic  r adiation  thr-ough  a perfect  vacuum.  It  is  a uni- 
vj'r-sal,  dimensional  constaiit  equal  to  2.907925  t 0.000004  X lo'^  cm  sec”*. 

Time  fixes  were  obtained  by  a straightforward  application  of  spherical  trigo- 
nometry. C’rrnsider  the  case  of  three  stations  (I,  J,  and  K)  that  receive  a signal 
as  shown  below; 


FIX 


Station  I receives  the  signal  first.  The  distances  between  stations  are  represented 
by  a,  b,  and  c measured  in  radians  and  between  the  stations  and  the  fix  by  X,  Y, 
and  /,  as  shown.  We  will  develop  formulas  to  find  the  angle  o and  tho  length  of 
side  X.  This  will  locate  the  fix  on  the  earth.  From  the  differences  of  received 
times  convertftd  to  distance  using  speed  of  light  propagation,  we  can  define: 
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A = X - Y 
B X - Z 

thu.s 

Y = X - A 
Z X - B 

F'roni  the  law  of  cosines  we  know  that: 

cos  Y = cos  a cos  X + sin  a sin  X cos  o 

cos  Z = cos  c cos  X * sin  c sin  X cos  3 

Substituting  in  3 the  values  of  Y and  Z from  2 gives: 

cos  (X  - A)  = cos  a cos  X ♦ sin  a sin  X cos  o 

cos  (X  - B)  = cos  c cos  X sin  c sin  X cos 

Expanding  the  left  side  of  4 gives  us: 

cos  X cos  A > sin  X sin  A cos  a cos  X + sin  a sin  X cos  o 

cos  X cos  B + sin  X sin  B = cos  c cos  X + sin  c sin  X cos  P 

Simplifying,  we  get: 

cos  X (cos  A - cos  a)  + sin  X (sin  A - sin  a cos  <>)  0 

cos  X (cos  B - cos  c)  + sin  X (sin  B - sin  c cos  P)  = 0 


or: 


tan  X 


cos  a - cos  A 
sin  A - sin  a cos  o 


and: 


tan  X = 


cos  c - cos  B 
sin  B - sin  c cos  p 


clearly: 

cos  a - cos  A cos  c - cos  B 

sin  A - sin  a cos  o sin  B - sin  c cos  |1 

To  simplify  let 

B = cos  c - coo  B 
and 

R cos  a - cos  A 
Substituting  9 into  8 we  get: 

H (sin  B - sin  c cos  P)  P (sin  A - sin  a cos  o) 

Also  from  the  figure  it  can  be  seen  that  o and  P differ  by  a constant  which 
we  will  call  C(JN  which  depends  upon  the  configuration  of  stations  used. 
Thus  we  may  write: 


82 


k 


11.  /?  ■ « - CON 
Thus: 

12.  R[sin  H - -lin  c (cos  o cos  (CON)  + sin  a sin  (CON))|  “ P sin  A - sin  a 

cos  a 

or: 

13.  sin  o (R  sin  c sin  (CON))  ♦ cos  o (P  sin  a - R sin  c cos  (CON) 

(P  sin  A - R sin  R) 

To  further  simplify  let; 

14.  S R sin  c sin  (CON) 

T P sin  a - R sin  c cos  (CON)  and 

Q P sin  A - R sin  R 

Formula  13  then  becomes 

15.  S sin  o + T cos  a Q 

Rewriting  15  in  terms  of  cos  o we  have: 

' 2 

l(i.  1 - cos  o - C}  - T cos  Q 

This  is  now  squared  to  eliminate  the  radical; 

17.  S^d  - cos^  a)  - 2QT  cos  o + T^  cos“  o 

Collecting  terms: 

1«.  ccs^  a (T^  ^ S^)  - 2QT  cos  o i (Q^  - S^)  0 

This  is  a quadratic  equation  in  cos  o.  It  is  solved  using  the  quadratic 
formula  to  give  the  results: 

19.  cos  o 2QT  t \'4Q^T^  - 4(T^  ^ - S^) 

2(T“  ♦ S^) 

I'sing  the  values  of  cos  u obtained  in  in,  we  can  find  .X,  the  distance  from 
Station  I to  the  fix  by  using  formula  7. 

The  system  f)f  sensors  consistently  "fixed"  the  incoming  sferics  within  the 

34 

proper  gef)graphical  area.  In  fact,  fre()uent  tests  by  Railey'  confirmed  that  the 
sferics  were  consistently  fixed  within  100  km  of  their  true  source  location.  As 
described  in  the  discussion  of  the  data,  the  datum  points  (counts  of  sferics)  are 
valid  for  a large  area  Isiundcd  by  a 10-degrec  latitude  and  10-degree  longitude 
block. 
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\:i.  niKiMiin  I)  vltlt^^ 


The  threshold  arrays  are  system  intensity  values  over  a regularly  spaced  grid 
that  represent  the  minimum  field  strength,  expressed  in  decibels,  below  which  no 
signal  may  be  located  and  recorded.  Figure  A2  shows  a portion  of  a threshold 
array. 
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Figure  A2.  Sample  of  a Threshold  Array  lor  Part  of  the  Eastern 
Hemisphere.  Values  are  in  dB  as  defined  in  Section  3 


The  threshold  an  ays  i<  >•  t!..  csults  of  a model  developed  by  the  U.  S.  Air 
Foret-  wi.h  the  use  of  ..n  IBM  3 .l-tiS  computer.  Figure  A3  illustrates  the  process 
used  to  develop  the  arrays.  There  are  two  threshold  arrays  for  1972,  covering 
the  periods  January  to  July  and  August  to  December.  Two  arrays  arc  necessary 
due  to  a decrease  in  the  number  of  sensors  in  the  latter  period.  The  model  inputs 
arc  information  about  the  sferics  source  (that  is,  field  strength  and  location), 
information  about  the  sensors  (that  is,  location  and  average  individual  sensor 
thresholds,  an  attenuation  table,  and  the  number  of  sensors  required  to  determine 
a fix).  The  model  inputs  and  the  model  will  be  explained  in  greater  detail. 

\3. 1 Mixli-I  liiiMii' 

The  source  information  included  sferic  locations  and  field  strengths.  The 
si.urce  locations  are  regularly  sf)aced  grid  points  in  five  ilegree  increments  over 
the  Eastern  Hemisphere  (see  Figure  A2).  The  initial  sferic  source  field  strengths 
are  arbitrarily  chosen  values.  The  reason  for  this  will  become  clear  in  the  expla- 
nation of  the  model. 
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Figure  A3.  Illu-stration  of  the  I^roce.srf  Fsed  to  Develop  the  Threshold  Arrays 

The  sensor  information  is  straightforward  and  requires  little  amplifying  dis- 
cussion. The  individual  sensor  threshold  values  were  determined  by  adding  the 
operator-set  threshold  to  the  basic  machine  threshold,  which  is  a constant  20  dB 
above  a 1 mVm  ^ field  strength  reference.  The  operator-set  thresholds  varied 
with  thunderstorm  activity.  An  incident  sferic  activity  increased  to  a point  that 
the  sensor  approached  saturation,  the  operator  increased  his  sensor  threshold  in 
1-dn  increments  so  that  only  the  more  intense  sferlcs  were  recorded  by  his  sensor. 
The  resultant  threshold  values  were  recorded  and  averaged  for  each  sensor  over 
the  two  six -month  periods  (Jan  to  July  and  Aug  to  Dec)  of  1972. 

The  signal  loss  from  the  source  to  the  sensor  due  to  attenuation  was  deter- 
mined by  using  an  empirically-derived  table  that  was  based  on  measurements  con- 
ducted for  the  r.  .S.  Air  Force  by  the  Denver  Research  Institute.  The  actual  sen- 
sor equipment  was  used  in  the  measurements.  The  main  factor  that  affects  propa- 
gation at  VI, F is  the  nature  of  the  boundaries  that  form  the  wave  guide;  the  earth's 
surface  and  the  D-layer  of  the  ionosphere.  Table  entries,  as  depicted  in 
Figure  A4,  are  averages  of  measurements  that  determined  day  vs  night  and  land 
* vs  sea  propagation  results.  The  empirical  table  gives  signal  loss  due  to  attenua- 

tion, expressed  in  decibels,  for  each  one-thousand  kilometer  increment  of  propa- 
gation path.  Table  values  are  additive  along  the  propagation  path  from  the  source 
to  the  sensor.  .Simple  linear  interpolation  is  used  when  necessary. 
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Figure  A4.  Sample  Table  for  Determining  .Signal 
lx)S,s  Due  to  Attenuation 


The  final  model  input  is  the  number  of  sen.sor  stations  required  to  give  a fix. 
.At  least  three  stations  are  required  as  a lower  bound  to  determine  a fix.  Any 
number  of  stations  may  be  used,  however,  up  to  the  total  number  of  stations  in 
the  system.  A four-station  fix  was  used  during  the  1972  period. 


\ 1.2  1 11.-  \l...l.  l 

The  threshold  arrays  represent,  on  the  average,  the  aggregate  sensor  or  sys- 
tem sensitivity  to  sferics.  As  may  be  seen  from  Figure  A2,  each  array  entry 
represents  the  signal  sensitivity  value  for  a 5-dcgree,  latitude-longitude  block. 

The  method  by  which  a single  value  is  obtained  will  be  explained. 

Consider  for  ease  of  explanation  a hypothetical  four  sensor  system,  single 
source  location,  and  threshold  array  values  as  shown  in  Figure  A. 5.  We  assume 
further  that  a four-station  fix  is  required  to  locate  and  record  a sferic.  We  will 
discuss  the  method  of  obtaining  the  (>1.4  dll  value  valid  for  the  block  centered  at 
54N’  I03K.  We  assume  a spherical  earth  and  choose  an  arbitrary  signal  field 
strength,  in  decibels,  and  a source  location  at  the  center  of  a given  block.  The 
path  from  the  source  to  each  sensor  is  then  determined  using  spherical  trigonom- 
etry. The  signal  loss,  in  decibels,  is  then  calculated  for  each  path  using  the 
attenuation  table  previously  explained  (Figure  A4).  The  signal  loss  is  subtracted 
from  the  source  strength  to  give  the  strength  of  the  signal  upon  arrival  at  each 
sensor  site.  'I'his  signal  strength  is  compared  with  the  sensor  threshold  value  at 
each  site.  Since  a ” four -station  fix"  is  required,  the  signal  strength  upon  arrival 
at  each  of  the  four  sensor  locations  must  exceed  the  respective  sensor  thresholds 
to  be  recorded  by  the  system.  The  source  field  strength  is  systematically  increased 
in  0.  1 dl?  inciements  until  an  associated,  received  signal  strength  is  reached  at 
each  sensor  site  that  exceeds  each  sensor  threshold.  The  lowest  such  value  for 
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Figure  A5.  Illu.stration  of  the  Method  Fried  to  Determine  a Single  Threshold-array 
Value  (1)1.4  dB)  using  a Four  Station  Fix.  Threshold  values  are  in  decibels  as 
defined  in  .Section  3 


the  45N  IO.tF  grid  point  in  Figure  4 is  the  01.4  dll  value,  which  represents  the  sys- 
tem threshold  value.  Other  geid  point  values  are  determined  in  an  identical  man- 
ner. F igure  AO  is  a flow-chart  depiction  of  model  calculations. 

The  threshold  arrays  are  used  in  the  normalization  of  the  sferics  data.  It  is 
necessary  to  develop  meaningful  relative  values  for  the  counts  of  sferics  over  the 
grid  since  the  sensitivity  of  the  system  of  sensors  to  the  collection  of  sferics  is 
different  for  each  grid  point.  The  statistical  distributions  used  in  this  study  were 
assumed  to  follow  the  log-normal  law.  That  is,  the  peak  field  strength  of  all 
received  sferics  foliow  a normal  probability  distribution  about  the  mean  (also  the 
median)  with  values  expressed  in  decibels,  that  is,  on  a logarithmic  scale.  The 
probability  (Q)  of  occurrence  of  a sferic  greater  than  or  equal  to  a given  threshold 
value  (TjT  at  a grid  point  (i,  j)  may  be  expressed  according  to  the  log-normal  law 
as 
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where  o is  the  standard  deviation  of  received  sferics  in  decibels  and  T is  the 
threshold  intensity.  The  grid  point  values  of  the  threshold  arrays  are  lower 

limits  on  the  log-normal  probability  distribution  used  to  normalize  the  counts  of  / 

received  sferics.  Each  sferic  count  at  a grid  point  will  be  divided  by  the  appro- 
priate log-normal  probability  distribution  value  for  that  grid  point  to  determine 
an  estimate  of  the  total  sferics  that  occurred. 
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Figure  AO.  Flow  C;hart  Depiction  of  Model  Calculations  Used  to  Determine  the 
Threshold  Arrays 
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